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Introduction 


partial  differential  equation*  play  a  central  role  in  may  arena  of  pbyalcs, 
engineering,  and  applied  outbanatica,  Existence  and  unlqueneas  theorem*  dare  been 
proved,  and  general  properties  of  lolutloea  have  been  studied,  for  large  claeeee 
of  pebbles*  for  partial  differential  equation*,  However,  explicit  exact  solutions 
of  prcblss:,  cwpwrwwtetw  sav  for  engineering  applications,  can  be  obtained 
only  for  relatively  few  prOblsnaj  and  often  the  analytical  form  even  of  these 
solutions  is  too  complex  to  be  useful  for  practical  applications.  because  of  this, 
conalderable  effort  has  been  devoted  to  the  study  of  approacUata  solution  acUrads. 
These  fall  mlnly  into  two  categories  t  — rieal  nethoda  and  asymptotic  aethoda. 

Bsc  suss  of  the  stiaulus  provided  by  the  development  of  high  speed  digital 
cosgutera  numerical  analysis  ban  nsda  tremendous  strides  la  reoect  years,  and  for 
may  prdbl—  involving  r^ial  differential  equations,  mmerlosl  nethoda  are  ideally 
suited.  Itr  ra*  j«trpoev*,  however,  these  net  hods  ere  Umreotlcel  or  even  useless. 
This  ie  particularly  true  when  one  la  prlwrily  conoernad  with  such  queetlon*  a* 
the  functional  dependence  of  the  eolutlon  on  the  parameter*  end  the  data  «f  Uw 
problem. 

Asymptotic  ar.Scds  have  been  developed  for  sane  type#  of  problems,  in 
particular  for  certain  ry Otises  Involving  e  parameter,  guoh  methods  provide  one 
or  more  terms  of  the  aaymtotte  expansion  (any  for  large  values  of  the  pareneter) 
of  the  eolntlon  of  the  problem.  Tbey  are  applicable  to  many  pent! ana  for  which 
enact  solution*  are  not  available,  end  even  for  problene  which  have  been  solved 
vssMTtly  it  often  happens  that  only  the  asymptotic  expansion  of  the  solution  is 
sufficiently  ilml*  to  he  useful  m  prectlocl  applications.  furthermore  it  is 
invariably  true  that  the  nethodc  which  yield  the  asycqtotlc  cspesnlcn  directly 
ur*  very  much  stapler  then  the  procedure  which  involves  first  finding  the  onset 
eolutlon,  thso  it*  aapsgtctic  aspanslon. 

This  report  la  davotad  to  Urn  study  of  a  eortcic  sicca  of  acpngtoelo 
method*  far  linear  partial  differential  equation*,  a  acmtrel  feature  of  such 
set  had*  la  the  notion  of  "roye"  vhlch  ere  curves  or  siretgh-  lines.  the 
rays  ora  of  fundamental  t^ertanee  beoauee  all  of  the  fttnstiona  which  soke  <p  the 
artcu  term*  of  the  upftotlr  etpaneloa  ana  be  show*  to  ectlafy  ordinary  tt.*.Vwdl 
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equoticns  along  these  curves.  Thus,  in  a  sense  the  net  hod  Is  one  vtt  ch  reduces 
partial  differential  equations  to  ordinary  differential  equations.  Often  vhe 
latter  can  he  solved  to  yield  explicitly  the  desired  asycptccic  expansions.  In 
seas  cases,  however,  the  ordinary  differential  equations  cannot  be  solved  explicitly. 
This  is  a  lialtitiou  of  the  nethod  which  Is  often  overlooked. 

The  historical  development  of  our  subject  Is  already  suggested  by  the  term 
"ray"  which  Is  u  control  Idea  in  "gecretrieel  optics."  In  our  study  of  the 

i-:  1  * . .  •  -  -  - -  ~r~ - - 

of  on  optical  protlsn  is  Identical  to  the  leading  ten  of  tfc**  csyaptotle  expansion 
(foe-  large  frequencies)  of  the  solution  of  an  appropriate  problem  for  the  reduced 
wove  equation.  Thus  tho  asyeptotle  aetbexi  shows  that  gscastrlcol  optics  Is  a 
first  approxientioo  to  "wave  cities".  But  in  addition  to  this  important  insight, 
the  asthod  provides  further  terms  in  the  asysgtatic  expansion.  These  terms  ere 
of  course  particularly  iapartont  in  regions  (such  ts  "shadow  regions")  where  the 
geasetrleal  optics  term  is  in*.  The  presence  of  small  disturbances  in  geometrical 
shadow  regions  has  been  called  the  phenonenon  of  "diffraction".  Thus  we  ay  soy 
that  the  asymptotic  theory  *t  the  reduced  wave  equation  yields  not  only  the  olcsslool 
geasutrloel  optics,  but  alec  n  new  "geometrical  theory  of  diffraction."  We  shall 
not  sttespt  to  eunemrlse  the  history  of  either  the  elaealenl  or  aotern  theory  here. 
■ow*ver,the  rafereneee  to  this  report  provide  so  outline  of  easy  of  the  contributions 
which  have  constituted  the  modern  development. 

Let  us  Indicate  vary  briefly  the  steps  involved  la  the  asymptotic  asthod. 

Mr  prrthlsms  which  eon  be  solved  exactly  ,e>ms(Limtlom  of  the  *  system  le  expansion 
of  the  solution  shows  that  it  consist?  *rf  a  sum  of  tense,  eeeh  of  which  is  on 
osysptotic  series  Involving  a  "phase  function"  end  ea  infinite  ssquenns  of  "nepUtude 
funetlcae".  tar  ecspltx  problems,  wo  thorn  fora  nssims  that  the  solution  is  also 
a  sum  of  such  scries.  By  inserting  amah  e  series  into  tbs  partial  differential 
equation  we  find  first  that  the  pbeet  function  satisfies  i  first  order  partial 
differential  equr'loo  which  eon  bo  solved  by  the  "msthod  of  eharacteristlcs. "  The 
ehereetorivtie  curves  ore  the  rsyt  which  we  here  nentlaned  eed  the  characteristic 
equations  will  bo  os l led  "ray  eqsatione".  thus  the  phase  function  eetiaflee  se 


j 


ordinary  differential  equation  alcdij  the  rays.  Wa  also  find  that  the  amplitude 
functions  satisfy  ordinary  differential  equations  along  the  rays.  In  order  to 
find  the  rays  and  the  phase  and  amplitude  functions  it  1$  necessary  to  specify 
Initial  conditions  for  ell  of  these  ordinary  differential  equations.  In  sme 
cases  the  Initial  conditions  are  a  direct  eonsoquenee  of  the  data  of  the  problem, 
la  ethers  the  data  *r»  hhtnlned  from  a  "canonical  problem."  A  canonical  problem 
la  a  problem  with  the  same  local  features  aa  the  given  problem.  It  la, however, 
sufficiently  alaple  to  be  solved  exactly,  rne  required  initial  uuuUiiOu*  or  too 
given  problem  ora  obtained  by  examination  of  the  nsyopbotlc  expansion  of  tha 
solution  of  the  canonical  problem.  Our  use  of  the  worl  "single"  In  eoaneetlon  with 
canonical  problems  Is  perhaps  misleading.  Often  th.  solution  of  a  canonic- 1 
problem  is  am  ambitious  research  project.  Once  It  it  found,  however,  (and  expanded 
asy^totloaUy)  It  yields  the  uwceseary  Information  to  complete  the  asymptotic 
:c'-tion  of  a  great  way  pr  mb  Terns  which  cannot  be  eolved  ?Actly.  Thus  the  asyeptotic 
method  also  provides  a  w&ds  application  for  e»ct  solutions  of  "sls^le"  prcblmm, 
banes  on  uddUlatml  motivation  far  studying  them. 

It  is  by  now  clear  that  the  asymptotic  method  i wolves  several  improved 
eamagblome.  It  Is  therefore  reasonable  to  mat  whether  It  con  be  proved  that  it 
does  indeed  yield  the  ovysgAotlo  expansion  of  the  exact  solution  tl«.  given 
problem.  So  germ  el  proofs  of  this  (hot  have  yet  beeo  glvea.  Isverthelsss  there 
Is  abundant  evidence  of  the  validity  or  tn  method,  the  evidence  la  obtained 
either  by  comparison  of  the  results  of  the  asymptotic  asthod  with  the  oayaptotle 
ovpniMlen  of  exact  solutions  (where  such  solutions  are  aval  labia),  or  by  comparison 
with  ouoarlaal  and  experimental  methods. 

Chopper  A  provider  the  east  complete  illustration  of  our  asthod.  In  that 
chapter  we  harm  attempted  to  provide  a  unified  ousamry  of  the  exist  log  literature 
am  the  asymptotic  theory  of  the  reduced  wave  equation.  In  chapter  1  we  apply 
the  seam  methods  to  Ih  well's  equations.  Moat  of  the  rasters*  there  are  the  seas 
aa  la  chapter  A,  and  the  render  whose  primary  Interest  la  la  cade*  eUmdlag  the 
theory  oaa  emit  chapter  K.  We  have  Inc  luted  It  bassmes  we  view  the  asymptotic 
netted  a*  a  practical  procedure  far  solving  problems, and  there  Is  a  Isrp  demand 

for  the  solution  of  problems  lmalvl*  tbs  elmetramgamtlo  field,  purthwmore  the 
vmntor  «he±-ct*r  of  thl-  field  lhnhse*  complicatin'*  which  do  not  erlae  In  ti~ 
study  of  the  reduced  wave  syatlcs, 

Thr<*mhoul  the  report,  vectors  sre  denote'  *y  capitals. 
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A.  Asymptotic  Methods  fo»*  ths  Reduced  Wave  Equation 
Al.  Asymptotic  Solution  or  .he  Reduced  Hsvc  legation 

Let  us  consider  s  reel  or  coo^ls*  function  »(t,  X)  which  utlsfles 
the  us vs  eduction 


1 


9  *  17U  \t  c  °- 


(i) 


Here  the  reel  wnlum*  function  c(X)  Is  the  prcpMStlos  speed  st  the  point 
X.  We  shell  look  for  s  product  solution  of  (l)  of  the  tom  »  •  f (tMf ). 
If  we  insert  this  for*  Into  (l)  end  sepurete  vnrlable*,«t  ohtsln 


r  set  X  ■  In  (2)  on 
o 

side  hy  *«*?.  then  (2)  yields 


(2) 


■gw  we  set  X  ■  X_  In  (2)  sod  donate  the  (eonetsnt)  enlne  of  the  left-heed 
o 


-o?  -  «*(tj/#(t). 

on  euhetltutlnt  ())  into  (2)  we  ohtsln 

J 

A  us  -r  n  »  0, 
o80O 

•qwstlon  (t)  for  n  Is  os lied  the 
he  fslshalu  waitloB.  It  le 
e  eoneuet  -efereea*  speed  eQ.  To  tei 
refmetloo  e(|  •  tjett)  end  the  WOOhtV 

1  \  *••<*> 


(J) 


<0 


lade  it 


tf  t.te  defies  the 


A*  ♦  tV(X)  n>  0. 


to) 
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The  constant  m  la  called  the  angular  frequency  of  the  solution  beeeuee 
two  linearly -lnoependent  solutions  of  (>)  «r«  the  periodic  functions 
c(t)  -  e'ltK  end  e(t)  -  e1*.  With  then  we  cun  ton  the  two  linearly* 
independent  product  solutions  u(x)c"lMtsnd  u(x)eij  fc.  Since  the 
eonplex  conjugate  of  every  solution  of  (>)  is  elao  a  solution  of  (3), 

s  Mk  a 

It  follows  that  every  product  solution  of  (l)  of  the  fom  u(x)e  is 
the  ewa^ties  conjugate  ox  a  solution  u;  »;sc  fv» -t 

v(t,X)  -  u(Xje*l"\  (6) 

Therefore  it  wifficee  to  study  solutions  with  *ne#tlve  tine  factor"  of  the 
fora  (6).  If  a  real  solution  v  is  required,  the  real  part  of  (6)  la  e»«ch 
a  solution. 

We  ahnll  mw  consider  the  solution  of  (3)  for  large  vel»*»a  of  k. 

We  begin  with  the  observation  that  when  n(X)  la  constant,  (3)  adalte  the 
plana  wavs  solutions 

u(X,*).t(X) a1"**.  (7) 

■are  the  fetation  vector  X  la  a  ml  or  ccaplea  vector  of  length  Of^.u 
sal  the  snplltuda  t  (R)  la  *  r*e>  nr  coanlett  eonataat.  in  feet  it  follow* 
froa  the  taurier  Ulagrtl  theoren  that  every  solution  or  (3)  with  sonata  1  a 
la  r  atparpoaiuoa  of  pleas  wave  solutions  or  Use  tarn  (7).  Ha  aspotaotlal 
a1®*’*  1*  called  the  pknae  factor  of  the  eolutUa  and  m  dell  cell 
ak*X  the  not.  %  aaalogy  with  {7)  wa  stall  soak  adotloaa  of  (3)  of  tka 
tarn 

u(X)  .  ttX.kfc1****.  (8) 

I 

! 

■ton  Inserting  (8)  late  (»,  sat  entailing  the  phene  fhetor  »  ,  wo 

oUaia 


t 


a* 
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A1 

k2£(^«)2  -  n2]  z  ♦  21^.^  ♦  ikz&»  ♦  ■  o.  (9) 

To  solve  (9)  for  laree  values  of  k  *-e  assume  that  z(X,k)  can  be  expanded 
in  inverse  powers  of  k.  It  is  convenient  to  -rite  the  expansion  in  terns 
of  (ik)  in  the  form 

-  _  .  (10) 

a(X,k).y  nM(X)(lk)-a  -  y  zm(X)(ik)  r*,  sm  •  o  for  m  » -I,  -8,...,. 
m«o 

Ue  have  used  the  sign  of  asymptotic  cam  lit  y  in  (10)  to  indicate  that  th«* 

aeries  must  be  an  asymptotic  expansion  of  t  as  k  -*  •.  This  means  that  for 

much  n  >  0 

*  n 

*(X,k)  -  £  *B(X)(i*r“  ♦  o(k“°).  (H) 

*s© 

By  definition  the  order  symbol  denotes  a  term  tor  which  11a  kD|o(k  )|  •  0. 
He  will  assume  thst  the  expansions  of  and  da  are  obtained  by  termwise 
differentiation  of  (10).  Upon  inserting  (10)  into  (9)  me  obtain 

21  («)1"“  {CW5-  n\*l  4  C2Vr^m  4  *m*J  4  \~l]  ~  °*  (A8) 

from  (12)  it  fellows  that  the  coefficient  af  each  power  k 
must  be  eero.  tor  m  •  •!  we  obteln 

[W*  -  •  0  ,  (13) 

einee  ^  go  for  •  •  •»,  ...  .  "U ,  ee  me  aeetae,  t0  i  o,  (13)  laade  to 

tha  ^*g^iaa  fw  •* 

(V^2.n8(X). 


A1 


fbr  a  •  0,1,2,...,  the  vanishing  of  the  coefficients  implies 

*  Eo  **  *  0  (15) 

and 

*  **  «  '  »  b  -  1,2,  ...  (16) 

Those  eqyetlono  ere  celled  the  transport  equations.  We  uill  see  that 
'■Q  can  be  obtained  by  solving  (15)  and  the  other  i*  can  be  deterBined 
successively  free  (16). 


t 
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A. 2.  Phase.  vave-frorte  srrt  r,.v». 

The  eiconal  equation  (1,14)  is  a  first  order  non-linear  partial 

differential  equation  for  a(x).  He  could  obtain  required  solutions  of  \ 

j 

(1.14)  by  applying  the  general  theory  of  firat  order  partial  differential  > 

equation**.  However, the  special  fora  of  (1.14)  enables  ue  to  tahe  a  simplified  ' 

(though  equivalent)  approach,  end  avoid  ten*  of  the  ccnpli  cations  of  the 

general  theory. 

tha  cur  face*  of  eoottant  phaae,  defined  by  •(*)  ■  ccnetant,  are 
called  wavefront* .  the  curves  orthogonal  to  them  can  bo  used  to  solve 

(1.14)  for  e(X) .  These  curves  are  called  ra>—  (In  tbs  general  theory 
they  are  called  the  characteristic  curves. )  The  equation  of  s  ray  nay 
be  written  in  terns  of  a  parameter  s  in  the  fora 

X  •  (i^.Hj.Xj)  •  K*).  (1) 

The  condition  of  orthogonality  is 


•  Xi^l  4  -  1,4,3.  (*) 

Mere  X(X)  is  sn  arbitrary  proportionality  footer.  Open  dividing  (8)  by 
X  sad  differentiating  with  respect  to  e  we  cfetalr 


• i 

40  *j  tr 

rV.=l’l&,v» 

S-mlr’U 

4  •  1,8,3.  (3) 

*W  (3) 

d.i4)  n*u 

l 

X 

-  (fS1)- 

«,(0’ 

(4) 

tr 

1 

s 

(8)  and  (1.14)  «iws 

*8ee  b),  Chapter  8. 


(5) 


earn**** 
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Equation*  (4)  are  a  system  of  three  accond- order  ordinary  differential 
equation*  for  t.he  rny*  X(n),  and  (?)  determine*  the  variation  of  the 
yeraneter  a  along  1  ray,  onee  X  haa  been  chosen.  «e  call  thaae  equations 
the  ffv  zt  la  to  be  noted  that  a  doea  not  occur  in  the*.  Hence 

the  reys  are  determined  solely  by  nW.once  Initial  values  for  (4)  are 
specified.  Of  all  the  rays,  &  two  porcssctsr  farily  are  orthogonal  to  the 

rtf  a  n)ir  n+  fnnrt.^  nw  *. 

If  ve  choose  X  *  n*^  the  ray  equations  take  the  fora 


Froa  (7)  we  we  that  for  thi*  choice  of  X,  w  ia  Juet  err  length  alone  the 
ruy.  If  w»  choose  X  -  1  and  denote  c  by  •<  the  ray  equation*  taka  the 
alapla  fora 


7r»  (7)  ard  (9)  it  it  alter  that  If  e  denotes  are-l'^th 


ft" 


\A  T* 


(io) 


eolve  the  *1  coral  equation  (1.14}  :  w  a  we  note  that  (1.14) 
and  (8)  yield,  fur  the  derivative  of  a  along  a  xey,  the  reaalt 


i 

f 

i 

j 

i 


.8. 


A3 


Is  «DC(«)].  **  %  -A  (*=)“->«*. 


do 


<n> 


Upoc  integrating  (11)  with  raapect  to  o  we  obtain 

•Dt(o)3  -  «Oc(o0)]  ♦  /  ADt(«,)]n2pt(o‘ado* 

®0 

When  X  «  n*A,  odm***#  arc-lenrth  and  (12)  baccnea 

i  f  -I \  «»<«• 

•v"  V 

Hera  ve  Hre  written  e(»)  for  •[*(<;)]  and  and  e  alallar  notation  far  n. 
SlaUarly  when  X  •  1,  (12)  becoaea 


(12) 


(13) 


a(f)  »  «(t0)  * 


£  n^T'Jdt'. 


OO 


(Ij)  and  (it)  provide  .iapla  fomilaa  for  the  value  of  ■  at  any  podafc  » 
a  ray  la  tana  of  tha  value  at  a  si  wan  fd«t. 


AV,  gnlwttan  of  tha  trapaoort  awU«Mjy  Vtt 

Za  tha  preceding  oaetton  tha  ray*  war*  uaad  to  cfetaia  tha 

aolutioa  a 00  of  tha  alcowl  aquation  (lUk).  nay  nan  alae  ha  naa*  to  aolva 

tha  tme*ort  equation.  (l.»)  «"d  (1.H).  «•  «**  “*•  “at  la 

HOfortlnnl  to  tha  dlraoUooal  darlvatlva  of  In  the  dlnetleo  of  %# 
wUah  la  Jut  tha  ray  direction.  In  fact  trm  (W)  «  obtain 

•  i  SI  *  *  X  fc  ‘nO^'O’  (l) 

fhu  we  •  .a  that  tha  traaaport  aquattona  (1.13)#  (I.U)  ara,  la  fhot# 
flrat  or«ar  ordinary  diffaraotlal  a«Utiaaa  aloof  tha  *aye,  MiufM 
wrlttoaaa 


(Jr  •■»*•. 0, 

(#) 

|  <^j|  ♦  a^hi  •  »  •  !#*»•••• 

0) 

w 
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We  will  first  obtain  the  solution  *Q  of  the  haeogeneous 
equation  (2)  and  then  um  it  to  obtain  the  solution  of  the  inhcaagenaous 
equation  (3)  by  standard  eethods.  Actually,  the  solution  of  (2)  la  aoat  ! 

easily  Obtained  by  returning  to  the  fora  (1.13)  and  noting  that  that  ! 

equation  lap  lie*  i 

| 

?•(»*  V»)  -  a  A**?1*  ♦  y*)  -  O*  00  j 

Given  a  ray,  we  nan  consider  a  region  R  of  X -space  bounded  by  a  tube 
of  rays  containing  the  given  ray,  and  two  segnenta  of  wave-fronts  W(oQ) 
end  v(o)  at  the  points  cQ  and  c  of  the  given  ray  (nee  figure). 


Thus  We  la  parallel  to  the  aides  of  the  tube  and  nomal  to  its  ends. 

Wb  new  apply  Uauas  theorwe  to  the  region  R.  gy  virtue  of  (k)  we 

obtain 

0uW  *{'i'*)am If '%*'**' If  •o*’**- 

*  WW  vTfo> 

■are  I  la  a  unit  vector  orthogonal  to  the  wave -fronts,  Bowsvar  froa 
(1.1k)  wa  aae  that  Pel  •  a.  tharafora  tqr  shrinking  the  Mbs  of  rays  to 

the  given  ray  wo  obtain 

•J(e)n(e)da(e)  .  aj(*0)n(er)aa(e0).  (6) 

tat  aa  aw*  choose  aa  arbitrary  point  ax  on  the  ray  and  eel 

<M  •  SS![)  (t) 


I 
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A3 


£(«)  la  called  the  eaaanalcn  ratio  alnea  It  neaaurea  tha  expend  co  of  a 
tube  of  raya.  It  la  Juat  tba  Jacobian  of  tba  Bopping  by  raya  of  V(o) 
on  W(o1).  Trtm  { C )  and  (?)  va  new  obtain  tba  aolntlon  of  (2)  la  tlu  tarn 


ftm  (8)  n  aaa  that  iQ(o)  varies  Inversely  aa  tba  square  root  of  n( 
along  a  ray,  ao  that  vban  (  dlnlniahaa  sQ  lncraaaaa.  thus  convergtnoe 
of  tba  *ty«  tanda  to  laeraaaa  a^  and  divergence  of  tbaa  taada  to  daeraaaa 
it  An  pbyalcal  Interpretation  la  perhaps  non  dearly  aaon  in  (6) 
ilcta  atataa  that  tba  energy  flux  a^da  la  ooaatant  along  aa  luflnltealaal 
tuba  of  raya. 

Za  order  to  obtain  tbo  aolutlon  of  tba  lnbcnoganaoua  agatln 
(3)  no  Intwodaea  tba  aolntlon 


of  tba  henogansoua  agnation  and  note  that  r(aQ)  *  1.  then  by  tba  netted 
of  "variation  of  pananters*  va  look  far  a  function  w(r)  aaab  that 

t^(<*)  ■  v(t)r(e).  (10) 

Zf  va  differentiate  (10)  vltb  regpeot  to  «,  lnaart  In  (3),  and  note  that 
r  satirise  (8),  va  Obtain 

S’-fe^Vl*  <U) 

Zt  follova  that,  of  do  ao  arbitrary  additive  ooaatant,  v  la  given  by 

*<•>--)/  <“> 

•« 
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end  the  general  aolutlon  of  (3)  1* 

C  i 

•m(o)  -  e1r(o)-*r(i)w(o)  -  ^rfo)-  x^a '^n-l^°  (13) 


By  letting  9  ■  oQ  we  aae  that  ■  ^(oq)  ,  hence 


■  »  1,8,...  .  (Ik) 


I*  we  chooeo  X  to  be  n"1  then  e  denote*  are-lcngth  along  the  ray  and 
X(o')  auit  be  replaced  by  n'^d')  in  (lk).  If  ve  chooee  X  to  be  x  then 
(lk)  beconei 


The  eolation  v(t,X)  of  (1.1)  ropreaente  a  dloturtanco  In  b 
phyaicol  eedlua  which  ie  eharcetorlaed  by  the  propagation  apeed  e(X)  or 
the  index  of  refraction  n(X)  *  Cg/cfX).  The  nediun  will  be  called 
{■xnorerroui  if  theae  function*  eru  oonataat.  In  thla  caaa  ur  earlier 
r.aulta  aiaplify  eonalderably. 

fir at  we  aee,  frun  (2. A), that  the  raya  art  etralght  linaa,  and 
frca  (2.1::  that 

*(«)  •  *(«*0)  ♦  nfo-Cp).  (1) 

gaze  a  denote*  arc -length  along  a  ray.  If  e  la  matured  on  all  raya 
free  a  warn -front  a(X)  ■  *Q  than 

a(o)  •  »(0)  *  m  ■  eQ  ♦  iw.  (2) 


» 


Benoe  the  distance  free  the  wave-front  »(x)  ■  *Q  to  the  wave-front 
a(X)  ■  »1  la  Just  (sj-a^/a.  But  this  distance  la  the  earn  on  every  ray. 
Therefore  the  wave -fronts  for*  a  fsmUy  of  parallel  surfaces. 

The  expressions  (3.8)  and  (3.14)  for  aQ  and  s^  can  now  he 


considerably  simplified  by  appealing  to 


elamntary  facta  of  the 


differential  gecaatry  of  surfaces:  Let  be  a  regular  point  on  the 
surface  8,  sad  let  H  be  the  unit  normal  vector  to  8  at  1^.  Svcry 
plsne  through  P^  which  is  parallel  to  8  cuts  8  In  e  curve  called  a 
noraal  section.  Let  *t  denote  the  curvature  end  ft  •  *_1  the  radlua  of 
curvature  of  the  normal  section  at  the  point  P^.  Then  <  depends  on  the 
direction  of  the  plane.  Zt  can  bo  shewn  that  there  exist  two  orthogonal 


directions,  the 


at  Plf  far  which  *  Lue 


and  minis**  values.  Theaa  values  are  called  the 
and  will  be  denoted  by  sad  *2>  Their  product 

«■  Yj’mT 


is  called  the 


of  8  at  P^.  Let  us  now  taka  8  to  be 


the  wave *front  V(s^)  (sea  aectton  AJ)  and  let  P^  be  the  point  of  Inter¬ 
section  of  the  ray  R  of  interest  and  the  surface  V(e^).  Via  plana  of  the 
following  figure  Is  chosen  to  be  the  plana  which  cuts  *(Oj)  in  the 


normal  section 


radlua  of  uunrature  la  p^. 


-1>  AU  j 

i 

The  ray  R  intersects  the  (parallel)  ware-front  W(o)  at  a  point,  P.  i 

Without  lose  of  generality  ve  nay  measure  e  from  the  > rare -front  W(e1) .  j 

i 

Then  -  0  and  the  distance  front  tc  P  is  0.  Since  the  ware-fronts  am 

per  a  lie  1  the  plane  of  the  figure  cuts  V(o)  in  s  narnal  section  with  radius 

of  curvature  ♦  0,  uud  u>«  ylau*  Uu-ouoh  K  orthogonal  to  the  plane  of 

the  figers  cuts  “(-)  is;  »  r->v»!e!  with  radius  of  curvature  0-  ♦  0. 

"  £ 

Purtheraors  it  is  clear  that  p^  ♦  0  and  pg  ♦  0  are  the  principal  radii 
of  curvature  of  V(c)  st  P. 

ia  the  center  of  curvature  corresponding  to  p^.  At  that 
point  the  ray  »  and  ita  nsidrtwinc  ray  IL ,  an  inflnltaalaal  distance 
■way.  Intersect.  Morn  precisely  ^  is  a  point  on  an  envelops  of  tbs 
rorilj"  of  re;-*.  There  is  s  rMUtr  point.  ^  correcpondHyr  to  the  other 
principal  radius  of  eurvaturo  Pg,  and  the  two  points  and  ^  lie  on 
a  ♦wo-eSeeted  snvelneo  of  the  ray  fasti  ly.  This  surface,  C,  Is  ealled 
the  etietie  or  mm?  fiTflBI  01  the  ray  tartly,  and  the  rays  are  taasont 
to  it.  ttons^lass  the  caustic  decensrataa  to  a  curve  or  a  point.  In 
the  latter  cast  it  ia  called  a  focus.  The  fatly  of  rays  itself  (each 
nomal  to  V^) ,  hence  to  all  the  wave  fronts)  la  called  s  gflQBl 
I MMEaMdd  rays. 

lew  let  b*  the  angle  between  the  rays  1  and  1^.  This 
ancle  subtends  ares  on  V(e^)  and  V(<?)  whose  Isaethe  are  P^dfl^  and 
(Pj+e)^  respectively.  dialler ly  we  jay  consider  a  ray  Eg  In  the  plant 
nomal  to  the  plane  of  tbs  figure  which  ashes  an  audio  ddg  with  S. 

The  latter  aa*l«  mbteads  ares  oa  V(?^)  and  W(o)  whom  lai*the  are 
Pgddg  sod  (pj+c)d«s  respectively.  It  follows  at*  (see  esetlon  A3) 
that  the  expansloe  ratio  is  dlwo  by 

i 

I 

\ 


gwt  g(o)  1«  the  Gaussian  curvature  of  V(r)  at  P. 

Iquttion  (t)  cna'ulesus  to  write  the  equations  (3.8)  and  (3*1*) 
for  the  amplitude  function  2q  and  in  the  alajxle  foru 

a0(a)  -  80(f0> 


V>  ■  V’o)£pi^)^  -  kj  <Vi,”')lk’'' 

®0  ■  .  1,2,...  .  (6) 

0l  and  Cg  are  the  prl-;:pel  radii  of  curviture  of  tha  wave-front  nt  d  •  0. 

At  thla  point  it  night  be  weU  to  point  out  the  conraction 

between  our  subject  and  the  mibject  of  eeonotrlcal  entice.  far  It  la 

clear  that  neny  of  the  teraa  vu  have  introduced  auoh  oa  "roar" ,  "ware 

front",  "oouatlc",  "fooua",  have  boon  borrowed  frea  that  afejeet.  A 

eloatr  rasper loon  ahewa  that  gscaatrlcal  rq>tlcB  ecnalata  of  a  aet  of 

rule  a  for  tho  cenatructleo  of  a  function  to  represent  wave  phamarna. 

lka(Z)  tw\ 

thla  function  tuma  cut  to  be  Identical  to  the  function  e  *qW 
which  la  the  leading  tei*  or  our  nsynptcti's  wapnneion.  (ftir  naarks  hare 
apply  aa  well  to  tho  emeu  of  inbewontneous  aedla,  n(X)  i  oonat. ,  aa  to 
the  caao  of  hcaonontoua  asdli) .  Indeed  tha  eayuptotlc  theory  oxplalaa 
the  apparent  p aratat  that  two  quite  dlfforont  baale  thaorlea,  fto» 

nptica  and  wave  optica  (i.e.  the  wave  equation) ,  have  been  ueed 
sm-oasafully  to  dnacribc  the  earn  physical  pitencaana.  Of  eoarao,  tha 
•ore  reatrietad  theory,  gnaaetrleal  optica,  can  be  uvoetad  to  ha  valid 
only  o*  high  ff?<g*enoiv»  (large  V) .  the  leading  torn  of  our  expansion 
la  often  called,  appropriately,  the  — c— trlcal  QQtlet  tan-  ftw  higher 


teroe supply  corrections  to  go  cec  trice!  optics.  In  recent,  years  U«r  sub¬ 
ject  of  geone tries!  optics  has  been  successfully  generalised,  principally 
by  J.  B.  Keller  and  his  co-workers,  to  explain  the  phencaena  of  diffraction 
and  other  phenomena  not  accounted  for  by  the  classical  theory.  He  vill 
hare  nore  to  ssy  later  about  this  uccnctrlcal  theory  fljfftpotlcn. 

Let  us  also  onelly  exsalne  Uw  cuuosction  Between  uur  suLjaul 
and  the  exact  theory  of  the  reduced  ware  equation.  Xn  special  costa , 
l.a.  for  special  choices  of  the  function  n(X),  and  for  prohlans  larolrl.^ 
boundaries  with  spsclal  gecnatricel  features ,  pr^blsu  far  the  reduced 
ware  equation  can  be  solved  exactly  by  the  natfcw?  of  separation  of 
ear  labia  a.  Man  ouch  a!*ut  solutions  ore  axpondrd  MysytotlcaUy  fer 
lam  k,  the  expansions  are  found  to  agree  exactly  with  cha  e^onaliwa 
we  ora  constructing.  Boworsr  Uw  class  of  prcblsns  which  eon  ba  solrsd 
exactly  la  extrenaly  restricted  and  even  for  this  clone  or  prehls— 
the  upptotle  solution  can  be  Obtained  nek  acre  quickly  and  aaslly  by 
the  present  nsthods.  Of  oourae  the  following  question  rsnalaai  far 
a  wry  wide  class  of  prefel  «s  the  solution  la  known  to  wist  and  to  ha 
unique,  but  solicit  construction  of  lU  exact  solution  la  w*  poaWble. 

Hi  theta  easts  Is  our  "aamytotm  expansion*  ladtod  ths  oaysytoUe  *x- 
pandlen  ot  ths  exact  aolutlent  At  preaont  no  pntwl  proof*  eon  ha  glean 
to  sowar  this  question.  Nsrertheles*  eensldsrdtla  xgyfw  ««d 
nwfrUTii  with  axset  solutions  which  can  bs  flhtnlnad  prerldsd  ws  with 
wanfldeacc  that  the  tMwr  e*  the  qsestion  la  afflie**! ee. 
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A5.  Wevee 

Ltt  ua  ramrlze  our  reeulte  for  the  etyeptotic  eolation  u(x)  of  the  re¬ 
duced  were  equation  (l.5)»  Fron  (1.-)  end  (l>10)  we  hare 

m 

u(X)  ~  ^  as(X)  (ih)*B.  ( 

The  pneac  inaction  aU)  ie  «  wuiution  at  tnc  .icwwi  «*u»U*u,  irA  ^.-twaing  to 
(a.lj)  ie  giren,  at  the  point  X(a)  on  a  ray,  by 


etXfe)}  •  st*(«c)l  ♦  J°  ntX(e')]ae'. 


Here  •  denotee  *re-len*th  alone  •  •*»«  *»•  ,r*  d^’tnioed  by  the  ray  equation 

(2.6).  the  amplitude  ia  given  at  the  point  X(e)  on  a  ray  by 

i(e.)nf>c(e ))  ^  .  . 

..w.n  ■  uw*.»^ j  ‘ 

and  the  other  ^(X)  are  qiree  rceurelvely  by  (3.1b). 

Hwb  the  fiwtlow  a(X)  and  i^(X)  hart  been  date  wined,  the  eeriee  (l)  ia  aa 
aagnptatie  aolutloa  of  the  reduced  rare  equation,  hah  a  eolation  will  be  called 
a  rare.  It  frequently  haaena  that  nore  than  one  ray  aaeooiatad  with  a  earn  paaaaa 
thrown  a  firaa  point  X.  In  each  eeaeo  the  vnlue  of  the  rare  at  X  la  Jim  by  a 

am  of  wptoaaiane  of  the  fom  (1),  one  for  each  ray  pnaaiatthrou#  X.  V  * 

«*tm  no  ttuonh  a  point  the  value  of  the  w**-e  at  that  point  la  wro.  ttne* 
tha  reduced  «w  equation  la  linear,  the  mm  or  etc  au*tr  of  eolutioae  of  It  U 
alee  a  eolation.  We  *al!  ere  thet  the  eajwptatlc  aolutloa  a f  a  «ireo  pruhlfli 
for  the  reduced  were  equation  will,  in  central ,  «»•**  of  a  mm  of  mraa,  tp- 
yroprUtely  aelteted  to  aatlefy  the  dot*  of  the  pretten. 


In  order  to  determine  c  vsve  uniquely.  Initial  values  for  s  Bust  first  be 
prescribed.  These  Initial  values  determine  the  fanlly  of  rays  associated  with  the 
wave  and  can  he  used  in  (2)  to  determine  s  at  every  point  on  every  ray.  In  addition, 
initial  values  of  the  functions  *B(x)  suit  he  prescribed  at  one  point  on  each  my. 
Then  these  functions  arc  given  at  every  point  on  every  ray  by  (3)  end  (3.1t).  The 
Initial  values  for  a  and  the  *_  ere  determined  by  the  data  of  the  problea.  for  ex- 
•■pie  In  a  radiation  problea,  in  which  the  solution  Is  generated  by  a  source.  the 
•ouree  will  determine  which  rays  occur  and  what  the  i«'tlal  values  are  on  than.  Xa 
a  boundary  value  problea  the  lnbonogeaeous  boundary  data  will  determine  the  initial 
values. 

To  recapitulates  A  solution  consists  of  a  eua  rf  waves.  A  wave  is  uniquely 
determined  by  prescribing  Initial  values  for  s(l)  and  Initial  values  foe  the  s^(X) 
on  each  ray.  Thus  the  first  step  In  determining  •  wave  la  to  solve  the  initial 
value  problea  for  a.  Ibis  is  the  subject  of  the  next  section. 


A6.  The  initial  value  nroblaa  for  the  eioonal  eeuatloa. 

Zh  the  usual  traatutwl  of  the  Initial  value  problea  for  a  first  order  partial 
differential  aquation,  Initial  values  of  the  solution  art  preaorlbad  on  a  surfaaa. 

Xa  ear  treatment  of  the  eleooel  aquation  we  will  also  have  to  coaaidar  lower  il- 
saulmal  initial  asai folds  specifically  curves  and  points.  Ms  will  mrtiw  these 
initial  value  problem  in  the  order  of  t  screes  lag  dinsweioo  of  the  initial  mmlfnll, 
i.a.,  point,  ajn,  an (  sarfaee.  ter  non*  of  thaae  probine  in  eha  aolntlnn  a(D 
indqeily  let  lied  >a<  by  the  value  of  a  on  the  Initial  mat  fold,  lowness ,  the  sain* 

tun  which  n  require  fur  c**r  Nvutmies  of  a  wave  u  eaiqaalv  letsrmlaad  bp  the 
additional  condition  that  it  be  *  outgoing”  tra  tsw  initial  analfnli 

A  aalatioa  e(x)  of  the  eteoeal  aqaattea  will  be  said  to  be  antgot.at  with 
ripest,  to  a  aneifhli  n  if  et  N  the  »ml  derivative  of  a.  te  •  ■,  in  positive 
ter  evssy  cutsets  aoraal  *  to  H.  iriiut  ftesttthen  every  dinettes  free  N  in 
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nor=?.l ;  If  M  is  a  curve,  the  normal  directions  at  a  point  lie  in  the  plane 
orthogonal  to  the  curve;  and  if  M  ia  a  surface  there  ere  two  normal  directions 
at  every  point,  one  on  each  side  of  the  surface. 

The  outgoing  condition  1*  consistent  with  t.he  physical  picture  of  e  dis¬ 
turbance  spreading  out  from  a  source  (whether  that  source  be  a  primary  one  or,  *t» 
in  ti.-  ctib*r  or  iBiAecti.cn  ty  f  *— •- *.^0:':  1-y  nnd  points,  a 

secondary  one) .  Mathematically,  the  "outgoing  condition”  is  the  aaymptotle 
analogue  of  the  radlstlon  condition,  without  which  the  exact  solution  of  c  problem 
for  the  reduced  wave  equation  ia  not  uniquely  determined . 

For  the  initial  value  problem  with  a  point  initial  manifold  P,  we  require 
an  outgoing  solution  s(X)  ■  hich  satisfies  th*  condition 

■(P)  ■  »0  .  (1) 

Clearly  the  solution  is  obtained  by  finding  all  the  rays  that  emanate  from  P. 

Then,  on  each  rsy,  a(X)  Is  given  by 

•  e 
'  o 

When  the  source  ia  a  curve  C  we  may  deacribe  it  parametrically  by  the 


s[X(o) )  -  a0  *J*  n  'X(a*))do‘. 


(2) 


eqvatior  ■.  •  XQ(T\)  ►■here  1\  denotes  srcl>.*ngth  along  C.  Lat  the  prescribad  value 

of  a  on  C  he  «[X  ft)i«  *„(?).  Differentiating  this  equation  with  respect  to  H 
o  o 


yields 


4X. 


da. 


(3) 


dX. 


7‘*  w  ■  W  ‘ 

lat  u*  introduce  the  angle  eft)  between  Vs  and  the  unit  tanfent  vector  ^  to  the 
curve  C  at  the  point  fl.  Then  since  the  length  of  V*  ia  n 

£  •  « 


coa  eft) 


nfXo0|)l 


t 


Because  the  direction  of  la  the  my  direction,  lift)  is  Just  ‘he  angle  between  a 

ray  leaving  the  curve  C  ot  T)  and  the  tangent  to  C  at  T|.  The  rsya  are  those 

which  emanate  from  the  curve  C,  at  every  point  along  it,  making  the  angle  p  with 

tX  tangent  to  C  at  the  point,  p  is  given  by  (U).  Thus  the  initial  directions  of 

the  reys  emanating  from  each  point  T\  on  C  lie  on  a  cone,  the  tangent  to  C  et  T) 

being  the  axis  of  the  cone  and  b(T\)  being  the  semi-angle  of  the  cone.  Then,  on 

every  ray  X(o;T\)  lying  on  the  conoid  emanating  from  the  point  Xq(yi) ,  a  ia  given  by 

a[X(a,T\)]  -  aft))  ♦  f  ntx(ff', J  dc*  .  (5) 

0  Jo 

da 

In  the  special  cnee  in  which  “  0,  J1*)  shows  that  p(T|)  ■  */2  ao  the  cone  la  a 
plane  normal  to  C. 

When  the  initial  nanif-d  la  a  surface,  3,  we  may  write  its  equation 
peraaotrlcelly  as  X  -  It  11  convenient  to  chooe*  the  parameters  and 

Tig  to  he  ere  length  a  along  orthogonal  curvet  on  R.  let  the  prescribed  value  of  • 
on  3  be  aiX^,^))  -  >0ftU  *V*  tifferentUtion  of  this  aquation  with  respect  to 

\  •*  %  *“Mb 

ax  ax 

gt  e  point  P  on  3  ,  ^  .*j  y£  are  orthogonal  unit  vector*  lying  in  the  teagent 
pint  to  8  et  P.  Let  P .  denote  the  angle  between  sad  gp  .  tnsn  (6)  fields, 
at  tha  point  F# 


««  •  i  w 1  i  “ 1,s* 

si 


These  equations  detemine  two  direction*  at  P  on  eppoeit*  sides  of  8.  These  nr 
t*»  possible  directions  of  '?».  Thus  two  rays  onenato  fren  each  point  of  d  on 


opposite  tides  of  the  surface.  On  eech  of  the  two  rnys  ana  ns  1  tag  fin 

Um  point  or.  3,  a  ia  given  by 

stxfol^n  •  *0tW  *j[  nlX^'JVV1  ^  w 
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For  completeness  ve  olco  Motion  th*'  charactcrletic  Initial  value  problem  for 

the  •  I <*00*1  equation:  We  sssubs  that  the  initial  values  s  satisfy  the  "surface 

o 

eiconal  equation":  * 

(9s0)S  «  n2  on  S  .  (?) 

Then  if  ve  choose  as  surface  co-ordinate  curves  the  level  curves 
*e  -  cuaat.  end  the  orthogonal  omdient  nurvee ,  ve  say  introduce  the  surface  para- 
««»:■  Tj  -  a  nsd  ?v  r.  is  sene  n**»»*t*r  that  labels  the  gradient  curves, 
e.g. ,  arc-length  along  one  level  curve.  If  Tj^  and  Tig  denote  arc-length  jmreiMtere 
corresponding  to  Tj  end  Tg  then  (9)  iaplles 
is  it. 

"  vq m  n  (10) 

and  clearly 

It  follows  froa  (7)  that  ^  ■  0  srd  «?  *  s/2,  i.e.,  the  rays  are  tangent  to  8.  In 
this  cote  since  8  <e  everywhere  tangent  to  rays  (i.a.,  rbaractarlstlcs)  8  la  said 
to  be  a  characteristic  surface.  One  outgoing  tangent  ray  XfoiTj.tg)  eaanatas  froa 
•vary  point  XQ(T1,tg)  on  8.  On  this  ray  s  la  glvan  by 

s(X(«jv1,t2))  ■  t,  +J  nWe'jTj.TgJldff'  .  (12) 

9m  surface  gradient  curves  fj  *  const,  are  everyimare  tangent  to  raye  and  nay  be 
celled  surface  rave.  They  play  a  central  role  in  our  later  discussion  of  diffraction 
by  saoetfc  bodiaa. 

*7.  Radiation  froa  sources 

(hte  way  of  characterising  a  source  is  by  giving  tot  values  of  the  phase 
function  :(<),  as  well  ns  the  amplitude  coefficients  ^(X) ,  at  every  point  of  the 
s«u ret  sent fold.  Usually,  however,  the  source  is  characterised  in  soae  other  way 
•nd  then  the  valuta  of  •  and  ife  oust  be  derived  by  procedures  vhlch  we  will  discuss 

•Tbe  surface  gmllent  $  la  OrfliMd  in  section  17. 


is  is  dr, 

O  0  4. 

3tj  A^g 


•  o. 


(U) 


1 

5 


i 
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shortly.  Let  us  nov  suppose  th.it  these  values  ore  given  and  examine  the  construc¬ 
tion  of  the  resulting  wave. 

The  phose  function  s(X)  as  veil  as  the  rays  are  d  etc  mined  by  the  procedures 
of  the  preceding  section.  The  amplitude  coefficients  may  be  obtained  by  mean''  of 
(3.14)  which  we  rewrite  in  the  form 

V0'  -  w 

D  9  0 ;X >2 fill  a 

Here  o  denotes  arclength  and  r.^(X)  ■  0. 

If  the  source  manifold  is  a  surface  S,then  on  every  outgoing  ray  we  may 
measure  a  from  S  and  ^(o)  is  riven  by  (l)  with  0Q  replaced  by  0.  We  are  assuming 
that  *m(0)  is  &tvn. 

For  point  and  line  sourcts,  the  source  manifold  la  a  caustic  of  tho  resulting 
ray  ayatam  and  banco  the  formulas  for  the  functions  ^  become  infinite  at  the 
source .  In  these  caaea  the  source  values  of  tVe  may  be  characterised  by 
appropriate  limiting  conditions.  These  conditions  will  be  given  only  for  the  cate  m«u.’ 

For  a  point  sourer,  let  dO  bo  an  element  of  solid  angle  of  tho  starting 
directions  of  the  rays.  Than  fer  sufficiently  small  *c,  da(oo)  *•  <rdD.  If  we 
Introduce  this  expression  in  (3 )  and  let  0Q-»  0  we  obtain 


da(c.)n(o„) 

da(c)n(a) 


-  a 


i 


TtTOori  ~~*wv~ 


uw  | 


#.  t 

Ui 


We  have  omitted  the  subscript  "o".  In  (0), 

«  (0)  -  Urn  9  t  (O.  (3) 

o-»  0 
c 

•Ve  will  assume  that  for  a  point  rouses,  ~  (0)  is  given.  For  the  eau  of 
hwncs norms  smota,  n  lc  sonercr/:  sri  d*(o3  ■  o^dfl.  Thar.  (2)  becomes 

The  analogous  formulas  for  m«i,d,...  sre  more  compl tcalv.l  mU  will  not  be 

require!  in  Uu? 


=  *i£L  . 

0 


For  the  source  diotributed  on  the  curve  C  (see  section  A6) 
it  is  clear  froci  the  foil  owing  figure  that  for  sufficiently  snail  ao, 

fo(a0)  ~  iIT]  sin  g  eea®.  (5) 


y,  sin  0 

/y 


Here  00  is  an  eleaent  of  angle  between  two  rays  lying  on  the  cone  of  rays  which 
annate  froa  the  point  TJ  of  C.  If  we  introduce  (5)  in  (l)  end  let  «0-*  0, 


we  obtain 


•(«)  -  *o) 


‘wfc-jf 


s  (0)  •  lias  s  (oj  . 


Again  we  sssum  that  ~  (0)  is  given.  Tor  the  cast  of  hooogeneoua  aedls 
we  can  ate  froa  the  following  figure  that 

da(«)  ■  ad®»(p,  «•  o)d*  •  M$(i  ♦  2-)sin  0  di).  (8) 

i  ox 


Bart  q1  la  the  distance  between  the  two  caustic  points  on  0*  ray  enanatlng  fven 
the  point  n  on  C.  The  point  t\  itself  is  one  caustic  point.  The  ether  point  T 
say  (as  In  the  figure)  He  on  the  backward  extension  of  the  ray. 


V 


-23- 


A7 


If  we  introduce  (8)  in  (6)  that,  equation  become* 


s(<0  = 


®(l  +  £  ) 


2(0). 


(9) 


-p^  is  the  signed  distance  from  the  curve  to  the  other  caustic  along  the  ray  in 
the  ray  direction  (i.e.,  the  direction  of  increasing  a),  this  distance  can  be 
found  by  deriving  the  equation  of  the  caustic:  A  variable  point  T  on  the  cone 
of  rays  emanating  from  the  point  XQ(t))  satisfies  the  equation 

(Y-Xo).Xo  -  |MJ  c°s  0.  (10) 

Here  the  dot  denotes  differentiation  vith  respect  to  1),  the  arclength  parameter 
w u  u.  Llfferentlation  of  (10)  with  respect  to  1\  yields 


(Y-Xe)-X0  -  1  -  -|Y-Xe|*  Sin  f»  .  .  X0eos  . 

(11) 

1 

3 

ty  inserting  (10)  iu  (ll)  ve  obtain  .  j 

(Y-X  ) *XB  9 

(Y-X0)-X0-  1  -  0  tin  9  ■  ag-fr  °  -  eeeV 

(12) 

* 

| 

or 

i 

(Y-X0)*(X0t  6  tan  B  •  sin*  0  . 

(13) 

1 

Ve  now  introduce  the  unit  tengant  vector  T  -  XQ  and  the  unit  normal  vector 
M  *  pXfto  the  curve  C  at  the  point  XQ(I)) .  p  denotes  the  rediua  of  curvature  of 
the  curve  at  that  point.  Then,  from  (lj)  end  (10) ,  the  eauatie  ip  given  by  the 
two  equations 


(Y-Xj*(ltp  B  tan  t  T)  ■  p  sin2  B 

(Ik) 

CY-X0)*T  -  |Y-X0|  cos  B. 

(33) 

SlUiuatiag  T)  from  (lb)  and  (13)  would  yield  a  single  equation  for  the  fanatic 
surface. 

low  let  b  be  the  angle  between  the  ray  end  the  vector  I.  Than  If  Y  is  the 


caustic  point  on  the  ray 


-<;4- 

(Y-X0)*N  -  .Blcos  ?,  (Y-X0)‘T  -  -  Px  ooa  p 
end  (lh)  yields 

a.  .  .  -  . 


A7 

(16) 

(17) 


g  fi  »in  p  ♦  cos  d 

In  (17 )  the  direction  In  vhieh  arclength  increases  along  c  Is  1  Material  since 
both  sin  p  and  p  are  unchanged  vh«n  T)  la  replaced  by  -T). 

We  conclude  this  section  with  a  discussion  of  radiation  frta  a  characteristic 
surface  8  (see  the  last  paragraph  of  section  a6).  In  thle  case  the  rays  are 
tangent  to  S  and  8  la  a  caustic  of  the  ray  aystea.  The  following  figure  'Shows 
the  "side."  "top,"  and  "end"  views  of  an  lnflnlteslaal  tuba  of  rays  leaving  the 
surface  Si 


(Vow, 


(for  sufficiently  iaall  oq  the  rays  arc  approxlaataly  straight.)  Hera  Sg  (r^) 
aeaaures  the  angle  at  which  rays  i 
Ola 

p2mvr  • 


its  from  the  curve  i  •  const,  end 
o 


free  the  figure  we  sae  that  far  sufficiently  mall  ^ 

«*(•„)-  90(e2*0)  r 

Xf  we  lotroduw  (l8)  In  (l)  and  let  «0-*  0  we  obtain 
•  (•)-»  (0)  jTa  sin  y  Sf§)j  • 


(10) 


s  (0)  -  Ua  e  ^  a  (•). 

0“»  0 
O 


(19) 

(80) 
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For  Uic  ctibc  Of  iiO*uOa*i'it?C,,.i*i  lucuid 


da(ij)  =  a(p2  +  c)  sin  y  dS^dd 


and  (19)  becomes 


>)  -  [o(l  4  H_)1‘ 
I 


For  some  purposes  (e.g.,  where  S  is  a  plane)  (19)  is  inconvenient. 

T*  bp  by  ’.'.sir®  (1 )  and  (w).  We  dei'ine 


dSfO)  da(jn) 

dSR  '  ““  0 

O 


z(o)  =  1(0) 


Idt(o)  n(0 


da(c)  n(<r 


A8.  Isotropic  point  source 

As  a  simple,  but  important,  illustration  of  the  foregoing  theory 
we  consider  the  problem  of  an  isotropic  point  source  in  a  homogeneous  medium. 

We  first  characterize  the  sou*®**  by  prescribing  the  initial  values  of  the  phase 
a  and  the  amplitude  coefficient  a.  The  coiidition  of  isotropy  implies  that  the 
limit  1(0)  of  (7.3)  is  ths  same  ir.  oil  directions,  i.e.,  on  til  rays.  Then,  if  we 
denote  distance  from  the  source  point  P  hy  r,  (-T.U)  yields 

s(r)  ■■  ,  (l) 

r 

From  (6.8)  we  have 

air)  «■  s(F)  ♦  nr,  (2) 


and  from  (%1) 


eU(.(F)w) 


If  the  source  Is  charaetcrlssd  by  ths  inhoaogsnsous  equation 
’/u  +  h8n2u  «  -ft(x-P';  fa  -  const,)  l  ( 


and  the  radiation  condition,  it  is  veil  known  that  the  unique  solution  is  the 
free  space  Green's  function, 
iknr 

u  =  T7T  •  (5) 

Comparing  (3)  and  (5}  ve  see  that  for  *  source  characterised  by  (4)  we  should  set 
s(P)  *  0.  ?.(0)  ■  rr  •  (6) 

The  problem  (4)  is  trivial  if  n  is  constant,  for  the  exact  solution  is 
i;ivcn  by  (' )  and  the  asymptotic  solution  la  unnecessary.  However,  if  n(X)  is  a 
variable  index  of  refraction  ve  My  consider  the  non-trivlal  source  problem 

V*u  t  k2n2(x)u  -  -ft(X-P).  (7) 

If  we  ass  urns  that  s(P)  u.-.d  K(0)  are  determined  only  by  local  properties,  then 

these  numbers  sre  given  by  (6)  and  the  asymptotic  solution  Is  given  by  (5*1),  with 

the  phase  given  on  each  ray  emanating  frost  P  by  (6.2) : 

s(c)  •a{x(o)J  •  f  n(x(o')]de',  (8) 

oJ 

and  the  amplitude  coefficient  s  given  by  (7*2): 


''  "  h  [' 


dfl  uCP 
d'a('c)  me 


The  problem  ve  have  solved  here  Illustrates  a  general  feature  of 
c-\r  asymptotic  method.  The  solution  of  tho  problem  (7)  was  determined  by  our 
earlier  conelderatlons  except  for  the  values  of  s(P)  end  l(o)  on  each  ray 
emanating  from  ?.  These  values  were  determined  from  the  exact  solution  of  the 
simpler  canonical  problem  (4).  Ve  shall  frequently  make  use  of  a  canonical 
problem,  which  can  be  solved  exactly,  to  obtain  certain  undetermined  coefficients 
for  e  more  difficult  problem  which  has  tho  ease  local  properties. 
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Agj  Isotropic  line  source 

In  this  section  we  examine  the  problem  of  an  iootropie  r.ource,  uni¬ 
formly  distributed  on  an  Infinite  straight  line  in  a  homogeneous  medium.  Let  r 
be  the  cylindrical  co-ordinate  measuring  distance  from  the  line.  Since  sQ  is 
assumed  to  be  constant  on  the  source  line,  (6.4)  Implies  that  p  *  jt/2,  i.e., 
at  every  point  on  the  source  line,  rays  emanate  at  rtp;hi.  angles  to  the  line. 

From  (7 .if)  wc  see  that 

.  .1 

1_  .  B  eln  &  ■»  0  cos  6  (,  \ 

P1  si n£p 

•  -1 

fo.'  3  -  0  and  the  curvature  p  of  the  source  line  is  aero.  Since  the 
medium  in  homogeneous  and  the  source  is  assumed  Isotropic  and  uniformly  dis¬ 
tributed,  the  resulting  wave  must  be  a  function  of  r  alone,  hence  (7 .9)  becomes 

a(r)  «  r"^  l(o)  (2) 


and  (6.5)  yields 

e(r)  -  s0  +  nr.  (3) 

It  follows  from  (5.1)  that  the  wave  produced  by  the  given  source  lo 

u-  I,lk(vnr)  •  (!») 

y/T 

Let  us  now  compere  (4)  with  the  two-dimensional  fret  space  Orson's 


function  ^  (knr)  which  Is  the  solution  of  (8.5)  in  two  dimensions.  By 
employing  the  euymptotlc  expension  of  the  Henkel  function  we  find  that  for  knr  -* » 

,n 

A  L. _ 

(?) 


$£)  (knr)  ._J -  * 

0  2v/5Snr 


and  we  see  thau  {'<*)  and  (?)  agree  exactly  If  we  take 


»  “o’0. 


i 

j 

i 

l 


«(o) 


(6) 
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We  may  new  consider  the  problem  of  an  isotropic  line  source  In  an 
inhomogeneous  medium,  characterized  by  the  two-dimensional  analogue  of  (8.8).  In 
this  case  X  =  n(X)  =  nfx^Xg),  P  =  P1»P2  6  is  the  two-dimensional 

delta  function.  The  leading  term  r>*  the  asymptotic  solution  is  obtained  by  setting 

1  lT 

s  =0  and  !  (0)  =  a  =  .  ■  e  ‘  . 

0  0  2/sxm 

Then  from  (fi.a)  we  see  that  on  every  ray  X  -  X(o;P)  emanating  (at  right  angles) 
from  the  source  line 


and  fro*  (7.6)  ve  obtain 

Since  n  &•  independent  of  jy  ell  rapt  twain  in  planes  ^  •  const.  Bene# 
da(o)  ■  dw(u)iUj  -  dv(o)dTV> 

The  — aniag  of  «w(o)  i»  most  easily  saan  frm  tbs  following  ffenret 
.  I  i  s 


by  &nssrtii%  (9  )  la  (6  )  and  collecting  our  results  ve  obtain 

"W**^*'  ♦<}  [si^feor]*  •  (l0> 


« 
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A  10.  He  flection  frca  a  bwafeg 

Let  us  suppose  that  i  wee  of  the  fws  (9*1)  Is  ImMbH  on  i  twisiy 
surlhee  1|  i.e.^the  nji  associated  with  the  wewe  Intersect  B»  let  the  solution 
n  he  required  to  satisfy  the  lMSdanoe  boundary  condition 

^  ♦  lks(X)  u  >  0(  X  on  B.  (l) 

Here  |>  >  Mfc  denotes  the  derivative  of  u  alone  the  aaftmurd  noewal  I  to  >#  sad  s 
la  s  given  function  called  the  fadeimi  of  the  boundary.  KrumndisO 
(l)  rednoes  to  the  ^er****  *  .-ludary  coulltlone  u  ■  0  end  ■  0  respectively. 

Hi  as  suns  that  s  reflected  wave,  also  of  the  fan  (9.1)  ie  yroteeed.  »»  verify 
this  ascent!  on  we  shall  show  that  (l)  osa  he  satisfied  hy  the  mb  of  the  tool  dead 
wave  u1  end  sa  appropriate  refleeted  wave  nr  which  we  shall  cosetrent.  thus  we 
write  a  as 

(ik>-eeU-r  £  ■^(ih)’* .  (») 

hh  mw  insert  (t)  late  (l)  eel  (Mala 


*  •  «  ♦ 


«r~  .“•*  «l 


/(X)  •  s1(x),  X  an  1, 


W 


t 


than  upon  aqaatlag  to  uro  tha  eoofflelaata  of  aaoh  powr  of  k  la  (3)  m  obtoU 


*)  ♦  (^  ♦  t)  .  0>  X  «  ^  (3) 

^(|*>)*  £(|$-+  •>  ♦  ^  -Ojn>ljXon».  (6) 


■potion  (4)  jrovldoo  tbo  mtaw  of  >r  oa  9,  (3)  oaa  ko  nM  for  a*  oa  M, 
aad  than  (6)  fatondaoa  anoeaoalNly  tfat  i£  on  B,  Xt  la  elnar  that!  thaaa 
mlaaa  onfflaa  to  datawdaa  tha  roOaotoi  novo  «r.  tot  «o  aa*  inalii  tho 
ofthUnaoa. 

tot  X  •  1||>  ho  tho  par— trio  motion  of  tho  Mahqp  oaa 

K  tha  raaatta  of  ooatloo  6  oonU  ho  naat  la  awlooMio  with  (t)  to  — 
tha  roOnetot  yfeaao  ftnotloa  ha*  It  to  —  no— law>  to  mnoooe 

tahtaMQyi  PUTaaaotUtUo  of  (h)  «&th  voavoot  to  aad  yUUi 


v,3fc***3!»  *“*• 


•njjl  aij^  UHikMaaioXttia^MaiM 

th1  u*  %T  mao  tho  am  promotion  on  tho  tanpmt  flam,  thnrafhai  tha  |l« 
aa— ag  thU  a— aa  inJnUsa  oof  tho  nmol  K  to  I  at  X  ooataiaa  both  ^ 
m«  V,  Xa  Jfman.al—  hath  a1  aai  a*  aatlafy  tho  la— 1  op—  fe.lt: 
th1  aaf  V  hoaa  tha  a—  laagii,  afe).  aiaao  that*  t— mill  oa— no—  m 


fe.lt). 


irthaku 


a—  at— Vo1  oof  V— MW 


that  or  ho 


ON  Of 


-■=?!- 


A10 


tVt  the  reflected  ity  dinettes  ^sr  lies  la  the  plea*  cctsladag  the  <«»«*«■* 
rey  direction  ^  ud  the  norael  I  to  B,  end  that  the  eagle  tf  refleettoe 
tf  equals  tbe  eagle  of  lirlrture  tf"  (See  figure) . 


the  egaellty  of  sad  0^  felloes  free  the  feet  that  the  aenel  eoe—e—e  of  \ 

1fcl  end  V  hen  the  sen  length.  j 

the  ledtlel  dlreotlae  of  eeeh  reflected  ray  Is  getetadaed  by  the  lev 
*?  -eflectlen*  sad  the  tadttel  nine  of  er  le  glees  tjr  (h).  «e  eev  eat 
a  »  a1  -  flf,  sad  nets  that  ■  a  see  a  sad  see  a.  thee  (9)  ead 

(6)  yield*  far  the  iattlel  eahtee  of  ajj  ,* 

« 

■J  •  K86b*  <•-«!  <•> 

vl  vf  (9) 

<*HS8  -infe. <%*♦%*>• •»*.«-». 


these  tadtlal  eelaee  eaahle  aa  he  i  asst  net  the  reflected  we(  Mae  eerlfytrg 
oar  aa— Hi  that  a  eolation  of  the  fan  (»)  setisflee  (I). 


All 


-?2- 

flusi  satisfy  appropriate  contiruity  condition*.  In  such  canes  if  a  save  Is 

incident  upon  S  that  surface  acts  as  a  secondary  source  and  produces  net  only 

a  reflected  wave  hut  also  a  transmitted  ware  on  the  other  side  at  tbs  interface. 

In  addition  If  the  boundary  or  Interface  contains  edges  or  vertices,  they  will 

act  as  secondary  line  and  point  sources,  respectively,  producing  what  we  shall 

can  diffracted  waves.  /Cu.  ana  multiply  reXlcuUu,  traaaivteS,  and 

diffracted  waves  nust  be  Included  In  the  sun  of  waves  fonOng  the  aoynptotic 

solution.  In  subsequent  sections  we  shall  shew  how  to  calculate  these  waver, 

A  11,  Be  flection  by  a  parabolic  cylinder 

Be  for*  considering  transmitted  and  diffracted  waves,  we  shall  1  Host  rets 

tbs  results  of  the  preco*4-’?  section  by  considering  the  probls*  of  reflection 
-ita 

of  a  plane  wave  ,  incident  along  the  axis  of  a  parabolic  cylinder,  from  the 


focus slz«  property  of  pemboUai  the  reflected  rays  are  redial  Hass  which  would 
|ms  throujfc  the  focus  if  extended  backward,  therefore  the  reflected  wave 
fronts  are  the  circular  cylinders  r  •  f7T?  ■  const. ,  i.e.  the  reflected  wave 
is  a  cylindrical  j»T2. 

In  genrrcl,  for  eylindrleol  waves,  one  priaslpal  redlae  of  curvature  of 
the  wave  front  r  »  const.  It  infinity  end  the  at  *er  equal*  r.  If  we  take  the 
wave- front  w(e^3  of  section  t  to  be  »  •  0,  then  o  •  r,  e^wo,  o2«  •, 
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All 


awl  (4.5,6)  beeooe 


r-*3 


(l) 


•  “o(v0,(r) 

\{r, 0)  -  *B(r0*®)(-f)^  -  ”^5  (r  )1^2e«*|^1(r  ,  0)dr;  ami, 2,., 


v  ; 


hy  ualng  (l)  and  (2)  ve  find  by  lndurfcioc  that 


«^(r,0)  -  | 

>•  ‘ 

a 

(3) 

Xaaartlag  (3)  into  (a)  ylalda  tha  raeuraira  fonulaa  for 

V0)"53 

(M 

(3) 

[vH[,*<w]1/‘ 

<«> 

la  the  valua  of  a^  at  mb*  polat  r0<®)  on  tha  ray  0  •  eoaat. 
Hb>  •  gyjjjrioal  vara  a  •  r  ♦  (*#  •  waat. ).  thwa  (9*1)  aad  (j)  ylaM 

fcb 

(7) 

teturnl*  dow  to  oar  laflaotloa  jrdbla»,wa  writ*  tha  aqpMtlaa  of 
tha  parabola  of  fooal  lao*U  p  aa 


* 

•/ 

1 


I 
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r  •  r.(0)  -  vJ?-.-, 


2  0 


rar©-pMC  2 


(8) 


On  the  parabola,  the  Incident  field  la  a"1***  e'^o  co*  Therefore  by 
lneertlag  the  total  field,  Incident  pin*  reflected,  into  the  boundary  coodltloo 

a  —  v,  aalala 

-lhr  coe  0  lks(r  .0)  •  r  ^ 

.  0  *  £  E4fo(0>»  3^)  -  0.  (9) 


If  we  equate  to  aero  coefficient*  of  pcwera  of  k  la  (9)  we  find  that 


a(v0)  • 

-r  (0)  00a  of 

0 

(10) 

*o[*o<0>» 

oj  -  .1 

(n) 

».[ro(e)» 

oj  •  Oj  *>  u 

(12) 

Aron  (10)  and  (8)  we  aee  that  an  each  ray 

a  *  aOr^O)  ♦  r-r^  •  r-r0(l*  eoa  0)  ■  r-2p,  (13) 

heoae  la  (7)* 

•0  •  •>».  (11) 


H>-  uslr-o  (11)  in  «)  we  cfctala 
foo(0)  * 


(15) 
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All 


tame*  fran  (3), 

=>,«)  - 


1 

2r  * 


(15) 


Ik*  we  my  determine  the  f  (0)  froa  (*»)  end  (5),  using  (12)  end  (15). 

By  calculating  the  firet  few  f  we  flsi  that  they  haw*  the  fora 

jb 

*  J^a#e  l,2J+1*  (17) 

and  (17)  can  be  proved  by  Induction  to  hold  generally.  In  (17)  the  a^  are  constants 
which  aatisfy  the  recuralon  fOnuiaa 


Vw»  * 11  (lfl) 

•on  -  -  ^  V  (19) 


•00  ’  -*• 


(20) 


Aran  (l8-20)  the  a  can  be  detendned  suoneeeively* 

Collecting  our  results  we  have,  for  the  aeynytetlc 
reflected  wave. 


esyanaloa  0 t  the 


tvi 


ik(r-flp) 


(a) 


problsn  treated  asyentotloally  In  this  section  oaa  he  solved  exactly 
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A12 


by  separation  of  variables  in  parabolic  co-erdiastes,  When  the  exact  solution  la 
expanded  asymptotically  for  large  k  It  yields  precisely  (21).  If  the  exact  solution 
Is  e capered  with  the  first  few  tens  of  (21)  (through  (kp)"2),  the  nuKriuel 
agraassBt  is  foutxl  to  be  good  for 

kp  >  2.  (22) 

lbs  prcibXeB  vs  have  discussed  herS|  and  nuasrous  related  prObleaa  are  treated  in 
[jtffl.  9m  dcaaln  of  validity  (22)  of  the  leading  tens  of  the  expansion  la  typical 
of  aost  problems  where  comparisons  with  exact  solutions  have  been  aade. 


A  12.  Reflection  sad  trsnsnlsslon  at  an  interface 

la  this  section  we  shall  determine  the  secondary  waves  which  are  prodarid 
whan  a  wavs  u1  is  Incident  oa  cos  side  of  s  surface  8,  across  which  the  index 
of  refraction,  n(x)  nay  have  a  juap  discontinuity,  flush  a  surface  Is  nailed  an 
interface.  We  denote  by  n^(X)  sad  n^X)  tbs  Halting  values  of  a  as  8  is 
approached  froa  sides  1  and  2  respectively,  and  wa  require  the  solution  u  to 
satisfy  the  two  boundary  conditions 


a1  •  au8,  ■  b  ,  oa  8. 


(1) 


the 


derivative 


■tie  a  and  b  are  given  functions  oa  8  sad  ■  >*^» 

to  8, 

We  s  $  suae  that  oa  the  aids  of  8,  say  side  1,  froa  whloh  the  wave  u* 
le  Incident,  e  refleeted  wave  ur  le  produced!  end  an  the  other  side  a  tnuaauittad 
ul  is  produced.  To  verify  this  aaeuardlea  we  shall  An*  that  (1)  is 
Mtlsfled  by 


a1  ■  a1  ♦  ar,  a*  •  at, 


(8) 
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vhtrt  ur  tod  u  art  aypi  opi'lataly  cooatnicttd  wvii;  outgoing  fro t  S* 

Thus  we  aet 

(3) 

u1  - X ^ •; (u)-, «*- u)-( 

where  *  i1'  •  ■  0  for  ■  ■  »  We  Insert  (3)  into  (2)  end  (2) 

■  me 

Into  (l)f  end  derive  the  folloslnc  conditions. 

84(X)  -  sr(X)  •  st(X);  X  on  Sj  (4) 

^  ■  es*  |  X  on  Sj  (5) 


If  X  ■  X6(ni  T^)  is  the  perenetrle  equation  of  '.Jam  interfS-ie,  8,  then 
different  motion  of  (4)  yields 

(t) 

mm*  these  eqnetlone  laply  thet  V,  end  V  4ns  the  sees  teojeotlw  ®  the 
1  ,  -  |0«ne  to  8  et  the  yolab  X.  It  f ullovs  thet  theee  three  rectors  end  the 
alt  eoral  vector  I  ere  eo*leaeri 


J 

N 
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JU.2 


ffcco  the  clconal  equation,  vt  jots  that 

c*1)2.^,  (*r)2.««J,  (tfr-4-  (e) 

It  fallowa  bow  mo  (7)^8), and  the  outgoing  condition,  that  9»*  and  we 
directed  as  shorn  1b  the  figure  and  that 


<?  ■  <?  |  Bj.  ala  9  *  n^  ala  (®4, 

(9) 

He  Mt 

0(^  ■  /  *8) 

(10) 

and  tins 

ala  P  ■  -I  ala  a. 

H 

(11) 

— 

(10)  le  the  fheUlar  lev  at  reflection,  aad  (U)  la  the  g£  j 

haturalag  to  (6)  vs  note  that 

^-■a^eosaf^>«-B^eosa^^  •  -Bg  oos  9. 

(18) 

Ms  laser* 

k  (12)  into  ($)  aad  lakrotee  the  ntlo 

hag  see  9 

(13) 

*  “  aa^  see  a  * 

9m  M 

aad  (6)  take  the  fora 

M*  -  *! 


X  oo  8, 


4«<-4*^4bi 


X  ob  8. 


tad  these  equations  an  easily  solved  to  yield 


r  jxt  rxl  «  (l£) 

i  *  fcr  *i  *TK)g^ooe5  £  “t1  *  Xon8» 

Ext  frl  ajr  -1<X7) 

*  “#J»  *  •  ** 


•i-itrfr'i  ^*TT5r 


(ib)  and  (17)  are  valid  for  a  ■  0,1,2;  ...  ,  tut  In  aneh  onae  the  eeeoad  ten 
on  the  right  aid*  veal  shea  for  *  *«,* 

the  initial  valuta  of  ar  and  a1  on  each  reflected  and  traantttad  ray 
are  (Ivan  by  (H),  and  the  initial  values  of  the  coefficients  a*  and  ■*  are 
given  by  (15)  and  (17).  iheee initial  valuta  stable  as  to  eoBatrast  the 
raflaetad  aad  trananittod  vans,  thus  verifying  our  initial  asswpticn. 

Xn  the  foregoing  discussion  re  have  tacitly  asavaad  that  (9)  can  bt 
solved  to  obtain  the  angle  0 f  refraction,  8,  for  all  angles  of  Incidence 
a  <  e/8,  ibis  is  certainly  true  provided  Bl/n8  <  1,  ■crevar|fcr  *]/*|  >  1 
than  la  n  critical  angle  gjKlflMt  a0  for  vblob  sin  a  1*  Mr  this 
angle  at  incidence,  the  angle  of  refraction  la  8.  ■  s/2  and  the  oorreepondlng 


angle  of  innidance,  the  angle  of  refraction  i 
treaadtted  »ay  it  tangent  to  ttw  Later  foe. 


aad  the  oorreepondlag 
e  for  a  >  a  (11)  has 

O 


no  real  solution  8  aad 


rliar  dlsousslao  mat  be  aodiflad. 


1  eoaplioetlGoa,  efeicb  occur  ebon  Bl/ag  >  1  are  sasocistad  nth 
of  ISEMtJIflflSllSfr  May  vlll  be  discuss  id  further,  lie  a  special 


>  in  the  folloring  aaatloa. 

m  ■  O.  tht>  111  mi A 


■  ’  0,  the  factors  “  end  nay  be  called  t»£LSiSiiM  «*>  jgaSiasMlgP 

o-fflri'-n’r. 


A1 3 
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&  13.  ggtoeMw  ind  teMdarta  of  a  cylindrical  gw  at  a  plane  interface 

I*#fc  n(x,y,x)  ■  for  y  >  0  and  n(x,y,s)  •  for  y<D,  vbere  and 
«2  ar*  firm  conatasts.  Cto  the  interface,  jiOki  prescribe  the  boundary 
condition*  (12.1}  with  a  and  b  constant,  The  Incident  war*  la  aanand  to  ba 


axeitea  oy  a  Use  source., 


fR.„  ",  •„  . . .  e*  •  •  1*  ..a ... .  .«  ai  . 

*«m»  miss*  ««•*  *«'  VttS  Vi  ViSC 


aourca  yroduees  tb*  Incidence  field  u1  •  |  I^fta^r),  (Sa*  aactlon  A.9)  tha 
tnnldaat  «m  la  (Iran  by* 


char* 


a^a^V 


»J(r)(Ui)*“|  y  >  oj 


(1) 


Baeauaa  of  tb*  eye-try  with  re  apart  to  tba  y-azla  va  aaad  calculate  all  function* 
only  for  x  >  0. 

Fran  tb*  law  of  reflection  it  la  easily  aaao  that  tb*  jtu*  ar  of 
tb*  reflected  wi  la  ar  ■  a^r'vbare  r  ■  /*8«(jr%)^*  la  tba  dlataoe*  fna 
tb*  yd  at  (0,-ai  rayr  aaad  lag  tba  "langa"  of  tba  aaaraa,  fbarefore  tba  roflaotd 

(1)  and  (2)  are  obtained  froa  the  aeyeptotlc  eapanalor.  of  the  Henkel  functlc... 
for  *  •*  *  m. 
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nn  li  given  by 


ur^eitalr  V  *1  (Ik)'* 


S 


(ik)*"  j  y  >  0.  (3) 


j* 

Mima  ♦Ha  w#1  •<»♦  o4  ww» .  ^«m»4  m  •««  <%4 «m»i! a «  mil 4  •*  Jj»y  Va 

- - -  m 

obtained  from  (11,3)  -  (11,6)  by  replacing  r  by  r  and  k  by  kn^.  We  also 

I 

set  r  ahieca  because,  for  points  on  the  Interface,  r  -  h  sec  or.  Thus 

zrjr\  a)  -  \r*  ^  yoHr'r^-i ,  (4) 

y°>  -  5j[(^)2  fj-i,A-i +  t*0**  '1*  w 

1  (6) 

tm( a)  -  (h  mc  a)2  t*  (h  — c  a,  a)  -  ^  (h  eee  a)"Jf^(a),  ■  >  i) 

f  (a)  -  »*(h  MO  a, a)  (h  eee  a)^Z.  (7) 

oo  o 

The  pbaae  of  the  incident  were  at  the  point  (h  tan  ato)  ie  e*  •  njh 
fharefore  the  pfaaee  of  the  tianaadtteA  wave  at  a  dietaaoe  e  fraa  that  point 
(along  the  traaaadtted  ray)  ie  a*  •  i^h  too  a  ♦  Cj  e.  Oau  the  tranoitted 


wave  la  given  by 


ut  ^  #ik(n^h  eeo  a  ♦  nj«) 


L 


pareaetrio  equation  for  the  trwuadtted  ray  ie 

X  ■  (x,y)  ■  (h  tan  a,o)  ♦  e(ain  P,  -  eoe  9)  (9) 


A13 
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and  rrOB  the  Uv  of  refraction,  (12.U.), 

•1  °2  ®1 
sin  3  ■  M  sin  at  u  ■  5=  -  -» 

*L  e2 


(10) 


Since  the  media  are  bcoogeneous,  ve  may  use  (4.6)  to  determine  the 
functions  Because  the  entire  problem  is  Independent  of  i,  ve  nay  take  one 

m 

radius  of  curvature,  say  pg  to  he  infinite.  In  determine  «  suet  flul  Ui 
caustic  of  the  transmitted  rays,  l.e.  the  envelope  of  the  family  of  straight 
Lines  (9).  Using  3  as  a  parameter,  ve  denote  the  caustic  curve  by  X  ■  X(3i. 
pj,  vlU  he  the  distance  from  the  point  (h  tan  o,o)  to  the  caustic,  alone  the 
backward  extension  of  the  transmitted  ray.  Thus  by  setting  *  ■  In  (9)  ve 
obtain 

X(p)  -  (h  tan  a,  0)  -  p^(*ln  3,  -  coa  3),  (U) 

and  differentiation  with  respect  to  3  yields 


(12) 

•  (h  sec8  a  H  ,  O)  -  p1(cos  3,  sin  3)  •  (sin  3,  *  00s  3). 

Sines  the  vector  |j|  is  tsngsnt  to  the  aeustlo,  it  Is  perellel  to  ths  rey, 
hones  perpendicular  to  the  vector  (cos  3,  sin  3).  Xt  follow  that 

(13) 

0  •  ||  •  (coa  3,  da  3)  »  h  sec2  a  eoe  3  Jj  -  p^  . 

let  from  (lO), 

toe  3  H  «  u’1  cos  a,  (14) 


Ox  ■  i»  h  eec3  a  coo*  3  -  u'H  eee3  oK»»*-eln*  a).  ^ 


*13 


Oils  rolae  of  &1  is  to  fee  used  in  the  relieving  equation,  obtained  trm  (1.6 )« 

&  «  C*®) 

sj(«)  •  (1  ♦  jjT  ^(o)  -  (1  ♦  -~)^2  de’)|  ... 


Ai  usual,  w  take  ■  0. 

The  funotioos  aad  «*  are  ftren  fey  (16)  end  (t-7)  oner  the  Initial 
tBlues  s*(o)  aad  s'  (h  see  a, a}  «re  specified.  Uammx,  these  talues  ere  dees 
fey  <12.16)  aad  (12.17).  Thu* 

(17) 

s'(h  •«  a,a)  -  ^  «i(h  *ec  a)  ♦  (141  "#*]»**»« 

r  ^  ^  -i  W) 

4(o)  -;fo  *i(h “° a)  ^Bcre^eora  [  if1  -  ‘  . 


a  e 


Ares  (8),  (b),  (7)#  tad  (17)  «e  easily  obtain 

.  -X,  -fcffr' )•*/*.  ./.’W,  <  .  «  M.  ^ 


Ififl  jw  ('^1  •**  (i*)( 


•  tri 


j 

4 

'0 


-U- 


(20)  iM  (21)  %Lrt  s*  and  explicitly  as  function*  of  x  and  y.  *o  obtain 
»*(x,y)  frca  (22)  it  would  ba  nacesaary  to  obtain  «(x,y)  and  <*(x*y)  fro*  (9)  and 


Zs  ordar  to  fbeilltata  tha  computation  of 


j£] 

TF*j»0f 


wfeieb  will  bo 


naodod  shortly,  It  la  ecmaolant  to  si^lify  (22).  «a  first  nota,  from  (9) 


..  «  f%  a\  it  .a. 

UM  \*‘V  V«BV 


*  "  Aa<**“  *“  a)  "  Ska  ‘eoa*a  * 


ftm 1  (19)  and  (23)  wa  bow  Obtain 


- 1  -  fa  *  db  [•»  ♦  -w*  *•*’ «]  •  £-*K“‘S)  } 


wa  now  insert  (19)  and  (at)  in  (tt).  tbia  yield* 


"o-OTI) 


li/a  * 

)  -at 


coa  Sf»Hi(iV 


ib»art ' 

» 


frtifl  1V*. 

a^Qr<b(l»u**)tani3  <gj 


Iron  (9)  wa  note  that 
r  •  Ni  I,  »*htan«*»ai«i. 


and  differentiation  with 


ta  y  yia'da 


*aao,  wbaa  y  -  0  (lfl),  (*)  aid  (?7)  i*Hr 


V 
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*  ■  s/2,  sis  a  •  u,  a  >  a  *  sin*1*!,  taA  «•  — ,  * 

0  *i-i 


(26) 

b^TM-vri 


tan  a 


He  nay  m*  uee  (28)  and  (25)  to  find  ^2-j^.  Dm  e<n*jtation  le  greatly 
airlifted  toy  noting  that  eoe  0  ■  0  when  y  ■  0,  hence 


Oi 


2  eln  a 


nV2 


eoe 


Bl 


(29) 


5-iywO  ‘  =.(!-** j  eoe  a  |_  ^fx^d-n^Hen^j  *  I*"0 

- X  (%)*/* 

i  -04  «0  o 

a*1  a»r 

The  derivatives  ^2  and  can  be  obtained  directly  fran  (20)  end  (21). 

Than  a£(h  eee  a/x)  and  **(0)  mi  be  obtained  fron  (l?)  and  (10)  and  need  la 
(k-f)  and  (16)  to  find  and  »*  along  their  reepeetlwe  raye.  He  will  not 
carry  out  thle  oalnulatlon  bare.  The  function  aj[  bee  been  maenad  in  §8J. 

If  it  <  1  total  rellactloe  occur*  for  reye  Incident  at  angina  a 

greeter  then  the  crltleel  angle  <*#  •  ale'St,  The  treaaadtted  reye  eorreeyood 

angle#  of  Incidence  a  In  the  Interval  0  <  a  <  cover  the  entire  lener  half* 

tjact,  and  the  critically  treneadtted  ray,  for  which  0  •  m/9,  Uee  la  the 

later  fee*  y  •  o.  Mr  reye  liwtldanr  at  engine  a  >  the  oorreeyandlag 

reflected  raye  are  called  "totally  reflected  reye*,  and  the  ecaTaej—dli 

of  refract  Ice,  0,  le  ecnglea.  «eeee  ae  real  treanadttad  reye  originate  le  the 

of  total  reflection",  n  >  h  tea  O^,  of  the  interface.  If  any  vnve  la 

le  the  eyec#  below  thle  regloa  It  la  pro  eel  la  the  aaee  ragioa  ae  the 

itted  vave,  whoa#  leya  originate  in  the  "renloe  of  renter  reflection", 

0  <  x  <  h  t-a  .  We  will  so  that  a  wave  la  yrodMced  belar  the  region  of  total 
•  •  o 

Tide  wave  e*  1<  called  the  eve— act  wave  beniiee  lta 
tepidly  with  dletaeoe  frna  the  intarfhee.  tb 


to 


with  the 
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tha  incident  vmve  u  end  the  reflected  ware  u  ,  satisfies  the  boundary  conditions 
in  the  region  of  total  reflection,  but  the  additional  presence  of  the  transmitted 
«nve  u*  belov  th_s  portion  of  the  interface  requires  that  e  fifth  wave,  the 
diffracted  wave  ud,  be  produced  above  tha  region  of  total  reflection.  Then  ud 

t 

together  vith  u  aatlsfiea  the  boundary  conditions.  Typical  rays  assoolatad 


Tha  total  field  vlll  bj  given  by 


ul  ♦  ur  +  u*1  I  y  >  0|  (30) 

u-v^+u*  ,  y<Oj  (31) 

d  A 

but  u  and  u  are  aero  in  the  region  of  regular  reflection. 

be  have  already  seen  hoe  to  obtain  ul6  and  tha  regularly  reflected 
pert  of  u^,  In  the  etna  vay  u*  and  tha  totally  reflaotad  part  of  «r  satisfy 
(IT)  and  (18)  (with  t  replaced  by  e),  bu.  the  angle  of  refTeotloa  9  lc  eoglag 
therefore  the  deternl nation  of  u*  belov  the  interface  le  different  free  the 
dutanlnation  of  u*.  Macs  u*  Is  not  r  squired  for  the  ealanlatlon  of  a*  , 
that  function  "an  be  found  as  before,  t,e,  ,(81)  renaina  valid  la  tha  region  of 
total  reflection,  Mr  «f»0  u*  is  required  to  detendaa  a^.  Nwever.ee  vlll 
not  daterdne  u*  here,  be  only  reaark  tint  it  daoayn  exponentially  vith 
diet  ana  e  fro*  the  interface. 


Dm  wn  ua  i»  coopletely  determined  by  u*  end  the  boundary  conditions. 
To  oaleulate  ud  we  ret 

and  write  the  boundary  conditions  (12.1)  in  tto  Tom 

u4  •  au*,  ■  b  y*  ,  yO,  x  >  h  tan  O0,  (33) 

9y  ineertin*  (32)  in  (33),  in  the  usual  way,  we  obtain 

(3V) 

e*(x,0)  ■  •t(x,0)  •  n^h  eee  oq  ♦  n^x-b  tan  aQ),  x  >  h  tan  aQ| 


*f(ik)‘ 


t  ike*  »— •  t 

—  5* 


(ik)’ 


(32) 


•*(X|0)  ■  e«£(x,0)  ,  x  >  h  tan 


yO,  x  >  b  tan  «0. 


(3?) 

(36) 


lb  (36)  we  have  weed  the  fleet  that  •  ^bt«(0,l)  ■  0  what  yC,  x  >  b  tan  afl. 
Xhie  rollon  Uwadlataly  fNe  the  dirwetlon  of  the  orltiaally  trenealttod  ray. 


To  calculate  u  we  sust  first  solve  the  initial  value  yaroblan  for 

e4(x,y)  in  y  >  0  «ain<  the  initial  values  (3b)  on  the  surface  yd.  tmomUag 

as  in  eoction  6,  we  easily  find  that  all  of  the  diffracted  rays  era 

to  the  erltlaally  rwflaeted  ray  end  leave  the  Marfhee  at  the  aafla  a  to  the 

o 

Me  also  find  that 

a  (Jf) 

a  •  n^(hey)  seo  ag  ♦  n^e-tbey)  tan  oJ|  y  >  0  x  >(h<y)taa  t*0. 


atm*  the  nave  front*  of  u4  are  plane.,  n  la  called  a  general  plane  nave,  The 
adjective  "general"  1*  required  ‘because  the  aoqpditude  ie  not  constant,  TO 
determine  the  amplitude  coefficients  t4  ve  note  first  from  (29)  that  «*(x,Oj  ■  0 
became  the  factor  cos  P  vanishes  vhen  P  -  x/2,  It  than  foUovs  firaa  (39)  and  (4.6)  that 
t4(x,0)  ■  0,  hones 

*g(*»y)  •  0.  (38) 

Vjyja  (4,6)  ve  aae  that  s4  Is  constant  on  each  diffracted  ray,  because  the  radii 
of  curvature  px  and  p2  are  both  Infinite  for  general  plane  aavaa.  Hence  free  (39) 

(39) 

u(x.v)  .  «J(x-y  tan  0^,0)  ►  a*J(x-y  tan  ctQ,  0),  y  >  0,  x  >{hty)  tan  a,. 

Since  ve  have  not  computed  »*  ve  cannot  use  (39)  to  ctotain  aj.  ■waver,  If  ve  set 

*■0  in  (36)  and  note  that  •  *d.(0,l)  •  ^  cos  ve  obtain 
.  .  ^*(x,0) 

•  ^rs0  -*r  *  <“> 


fc*(*,0) 


end  (39)  and  (40)  yield 


4,  j  b  *0  <«V?  |y‘  .0 

a^(z«y)  •  ^  eoa  v  L*  ■  *-y  ten  a#* 


Finally,  ve  Insert  (9)  la  (*1).  the  result  ie 

•S'1  *  ,  >  0(  , 


V 


The  tens  •“*  (U1^  la  the  loading  '  era  of  «*  Tt  la  of  order  l/k 
et  to  the  lscideot  field,  end  of  the  awa  order  as  the  tern  uj,  ehteh 


%  ** 


AJ3 


I  A 

n  han  not  determined,  zj  becomes  infinite  on  the  diffracted  ray  which  coincides 
with  the  critically  reflected  ray, 

x  •  (h-ty)tan  aft.  (43) 

Therefor*  the  preaewt  arvK/tcrtlc  expansion  fall*  on  that  ray. 

We  note  that  the  diffracted  wave  la  coc^letely  determined  by  Its  Initial 
values,  all  of  which  are  given  by  (3?).  Bence  ell  the  quantities  in  (36)  ere 
already  determined  and  therefore  that  equation  must  be  an  identity. 
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A  lU  Diffraction  by  Edges  and  Vertices 

A  surface  or  curve  is  regular  at  a  point  If  it  can. 'be  represented 
by  functions  vhlch  have  derivatives  of  all  orders  In  a  neighborhood  of 
the  point.  An  edge  is  a  curve,  on  a  boundary  or  Interface,  vhlch  fanes 
a  locus  of  points  where  the  surface  le  not  regular.  A  vertex  le  an 
Isolated  point,  on  a  boundary  or  Interface,  where  the  surface  la  not 
regular,  or  an  isolated  point  on  an  edge  where  the  edge  Is  not  regular, 
teaples  of  edgee  and  vertices  aret  the  edges  and  vertices  of  a 
polyhedral  Interface  or  boundary,  the  vertex  of  a  conical  Interface  or 
boundary,  the  edges  of  an  spertare  In  a  thin  screen.  (in  the  last  oaea 
the  aereen  la  a  boundary  eurfaoe,  but  both  aides  of  the  screen  being 
connected  by  the  aperture  ,  foam  the  doaain  of  the  problse.  If  the 
aperture  edge  la  not  regular  it  contains  vertex  points. ) 

We  have  already  Motioned  diffracted  warco  and  diffracted  rare. 
w»  will  use  these  tana  to  include  all  waves  end  rays  not  predicted  by 
the  classical  theory  of  gecMtrioal  optics.  Whan  any  wave,  u1  la 
incident  upon  an  edge  or  vertex,  H,  we  Sanaa  that  N  acta  as  a  secondary 
source  eanlfold  producing  a  diffracted  wave 

«**•““*  jr  al<lhr«.  (1) 

if  aaalogy  with  our  results  fur  secondary  waves  produced  by  reflection 
and  t  tenant  as*  cn,  we  as  ease  that 

«4(X)  •  e4(X)  ,  X  on  M,  (2) 


where  s1  1b  the  phase  of  u1.  The  point  or  curve,  M,la  a  caustic  of 
the  diffracted  wave;  hence, as  we  nave  seen, the  functions  are  Infinite 
there.  However,  the  limit  introduced  In  section  A7, Is  finite.  We 
BBBume  that  la  proport  lorai  to  the  amplitude  s*  of  the  incident 
vcvm  ▼,  l.e. , 


^(X)  .  (d)  «i(x)|  X  on  H.  (3) 

The  proportionality  tnc„;r  (d)  will  be  called  a  diffraction  coefficient. 

It  le  analogous  to  the  reflection  coefficient  r  »  ^  444  >  which 

appears  In  (10.6), or  the  tiwtiitw  coefficient  «  J 4»  j » *hlch  appears 
in  (12.1?)  for  ■  ■  0.  In  general,  diffraction  coefficients,  unlike 
reflection  and  transmission  coefficients,  cannot  be  obtained  directly 
from  the  prescribed  boundary  conditions.  Instead  they  are  cfctalned  either 
fran  the  solution  of  oeaonlnal  problems  or  by  boundary  layer  methods^.  The 

A 

latter  eetbode  also  yield  the  *Ouee  Of  the  1^  for  a  >  0.  far  mm  purposes. 
It  sight  suffice  to  determine  (d)  eaperlMnully,  nut  this  baa  sot  yet 
been  attempted,  lu  general,  the  diffraction  coefficient  ihpes*  on  the 
looal  geoi.  trie  properties  of  M,  the  local  values  of  the  Index  of 
refraction,  the  directions  of  both  las  Ideal  end  diffracted  laye,  ami  the 
wave  nutter,  kj  and  It  antttsi  la  the  Unit  k  -*•. 

The  phase  s*(x)  and  the  rays  of  tbs  diffracted  save  are  obtained  by 


solving  the  initial  value  problem  for  the  el  coeel  agnation,  with  Initial 


tcIum  given  by  {?).  Sirre  this  tea  been  discussed  in  detail  in  Section 
A 6  we  need  only  aent^oo  the  consequences  of  the  special  font  (2)  of  the 
initial  values  when  M  ie  an  edge.  Let  us  first  asstas  that  the  index  of 
refraction  is  continuous  in  a  neighborhood  of  the  edge,  as  is  the  oase 
when  the  edge  lies  on  a  boundary  surface.  In  this  oase  it  follows 

tt  n~l  /el  **••* 

. “  “  ’  "  to 

eo,0(n)  -spftnjr  ft-  *1,  rH“00*aCn)* 

Im  d(n)  is  the  esui -angle  of  the  o«ne  of  diffracted  rays  as  lasting 
fra  tbs  point  XQCn)  of  the  edge  and  oc(n)  is  the  angle  between  the  incident 
ray  and  tbs  odgs  at  that  punt.  Since  both  angina  Ua  between  sero  and 
«  it  frillies  fron  (b)  that  they  are  equal.  Hub  we  have  obtained  the  mfei 
Jp-  yf  ^1“  eiffteetiaai  tbs  angle  of  diffraction  is  equal  to  tha  angle 
of  InoitaKo.  Incident  and  diffracted  rays  in  the  neighborhood  of  a 
typical  point  m  an  edge  are  illustrated  in  tbs  following  figure. 


If  lbs  edge  line  on  an  inter face  there  say  be  two  or 


la  the 


seah  wedge  bet  «t 


of  tbs  edge  with  values  of  a(X)  eoatlaeme 
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>n&  n4  and  n1  are  the  values  of  the  index  of  refraction  tn  the  regions 
the 

containing 'diffracted  and  incident  ray,  at  the  point  of  diffraction. 

Tor  vertices,  which  are  secondary  point  sources,  (2)  has  no  special 
consequences.  Diffracted  rays  eaenate  frca  the  vertex  in  all  outward 

iircftisr.c  in  thi  ics?!??  **  t***  Tnv»M<N* 

A 

Once  the  diffracted  rays  anil  U»  values  of  1  on  M  are  duteralned, 

■ 

the  s*(Xl  and  hence  the  diffracted  wave  u*  follow  laaedlately  frca  the 

u 

formulas  of  section  At. 


To  illustrate  tne  roregoing  theory  v«  uiioll  vwi»iu*i'  iose  prcblcr: 
vlth  edges  on  the  boundary  of  a  medium  vitfc  Index  of  refraction  ml. 

Wc  begin  vlth  the  case  in  vhlch  the  edge  is  a  straight  line  and  the  incident 
rays  all  lie  in  planes  normal  to  the  edge.  Then  the  diffracted  rays  are 
also  normal  to  the  edge  and  emanate  frcm  it  in  all  directions.  Thus  it 
suffices  to  consider  all  the  rays  in  one  plane  normal  to  Um  edge.  IT  i 
denotes  distance  frcm  the  edge,  then  the  phase  s4  of  the  diffracted  vave  is 
equal  to  i1  ♦  r,  vhere  a1  la  the  phase  of  the  incident  vave  et  the  edge. 

The  edge  lice  on  an  Incident  wave  front, hence  a1  is  constant  on  the  edge. 

Since  the  diffracted  vave  la  cylindrical,  s4(r)  is  given  by 

a4(r)  -  #(0)^2  •  (a)a1  r"l/8.  (1) 

Karel d)  denotes  a  diffraction  coefficient  and  a4  is  evaluated  at  the  edge. 

Thu*  the  leading  tens  of  the  diffracted  vave  le  given  by 

«,V**  (d),J  r-V*  (!) 

la  (?)»  u*  ■  elk*  s1  denotes  ths'gratttrleal  optica  term”  l.e.^ths  leading 
tara,  of  the  incident  vave,  evaluated  et  the  edge. 

fl*al 

let  us  onagers  our  result  (2)  vlth  H— rf eld* e  exact  solution 
far  diffraction  of  a  plant  vave  by  a  half-plane.  That  result  consists  of  the 
Incident  and  reflected  veves  of  gsonstrionl  optics  plus  a  third,  or  "‘diffracted* 
tens.  When  u«  third  tern  it  aapgtctlcaliy  expended  f w  »i«*«  of  ter 

it  e^eee  perfectly  with  (i ) ,  provided  that 

«Moet  cf  the  eeterlal  in  this  section  is  adapted  from  [ldQ. 


A15 


4r  .  , 

(A)  -  - S-tt 5 -  [sec  |(0-a)i  esc  «♦«)]•  (3) 

2(Z«k)  '  sin  0 

Here  0  Is  the  angle  of  incidence  (or  angle  of  diffraction)  vhich  is  s/2  <n 
the  case  ve  are  considering*  Kic  angles  between  the  incident  and  diffracted 
rays  and  the  normal  to  the  lialf -plane  are  a  and  fl  respectively.  They  arc 
illustrated  in  tne  following  figure.  [Thk  *eug«  is  *  nau-pisnc  vr.au  ;  - 


The  upper  sign  in  (3)  applies  when  the  boundary  condition  on  the  half-plane 
is  u  -  0,  while  the  lower  sign  applies  if  i*  1"  -  °- 

The  agreement  between  (2)  and  the  eaact  solution  of  the  canonical 
problem  (i.e.,the  ftsmawfeld  problem)  la  a  confirmation  of  our  theory  «d 
also  determines  the  edge  diffraction  coefficient  (d).  Wallar  agreement  occurs 
for  oblique  incidence  on  a  half-plane  uhro  (2)  ia  replaced  by  the  appropriate 
expression  and  the  denomimtor  tin  (  ia  included  ia  fl).  m  this  caac  0  end 
a  are  defined  as  above  after  first  projecting  the  rays  Into  the  plane  normal 
to  tie  edae.  In  caaa  the  half -plane  le  replaced  by  s  wedge  of  angle 

7  ■  (2-q)»  M 

caparison  of  (2),  and  its  modified  fore  for  #  /  s/2,  with  geamerfeld'e 
•set  solution  for  a  wedge, yields  apmement  when 

w  •  \sSkr»  E(“  *•' *(“ *•  “* ^  ” 

for  q  •  2,  the  wedge  becomes  a  half -plane  sod  (^)  reduces  to  (3). 
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We  ehall  now  apply  (2)  and  (3)  to  do l* mine  the  field  diffracted 
through  an  Infinitely  long  alit  of  width  2a  in  a  thin  screen.  For  simplicity 
we  shall  assume  that  the  incident  field  is  a  plane  wave  propagating  in  a 
direction  normal  to  the  edges  of  the  slit.  Then  we  can  confine  our  attention 
to  a  plane  normal  to  the  edgee.  In  thle  plane  let  the  screen  lie  on  the  y 
axis  of  a  rectangular  co-ordinate  system  with  the  edges  of  the  slit  at  x  »  0 
and  y  «  ♦  a.  Let  the  incident  field  be  the  plene  wave 

ui  _  eik(x  coe  a  •  y  ein  a)  (6) 

too  singly-diffracted  rays,  one  from  each  edge,  past  through  any  point  F. 
thus  the  leading  tern  of  the  singly-diffracted  field  at  P,  u^(P)  in  the 
«"■  of  two  terms, 

ikfrj-a  sin  a)+lr 

UJ(P)  i — C-«  i c,c 

lk(ra*s  sin  a)*ijj 

*  *  9^'yW - 8(  1  c*°  »( V^J *  (7) 

In  (7),  rx  sad  rg  denote  the  distances  tram  P  to  the  upper  sod  lower  edges. 


ef  these  ray*  begin,  at  one  edge  of  the  tilt,  la  dlftaeted  tram  the  other 


-‘>7- 


sdge,  and  then  passes  through  P.  To  find  tba  carrsbv*~ulng  wm  It 
1*  ue^esaary  to  treat  the  two  singly-diffracted  wave*  aasnetlng  froa 


1b  uecesaary  to  treat  the  two  singly- diffracted  wave*  eaanetlng  froa 
th*  tvo  edge*  as  new  incident  waves  on  the  opposite  edges  of  the  slit,  «nJ 
:.v  y  (?)  (?.)  vith  c*  -  «/?.  TV*  *<*ir»ifcat.ton  Is  etxeigfat- 

forward  when  the  boundary  condition  is  u  ■  Ot 

Proa  (7 )  v»  ■**«  that  at  edge  ( J )  the  leading  tens  of  the  singly  * 
diffracted  wave  sweating  fro*  edge  (f )  is  given  by 


u(ij  -  -  — 


ika(2>  siitiO+l  n/k 

*7^ - 


[sec  ccc  5(5  -  a)]  •  (8) 


Bare  we  have  used  the  upper  sign  in  each  tens  of  (7)  corresponding  to  the 
boundary  condition,  u  e  0,  and  have  chosen  the  appropriate  valuee  of  r^ 
and  Oj.  (8)  is  easily  simplified  to  yield 

.Uate-(-l)J  tin  oQtu/b  .  (9) 

u(4)  ~  -  **  fl9* 


Ttm  (8),  tbs  leading  tew  of  the  doubly-dl ffwcted  field  et  f  la  gis 


u!(l I  -  51  (d)  «(j)r7  C^J, 
J  J 


i*A  .  ,  i a/t  . 

<4)  *  *  5^?/*  ***  4  *  *  “fejT/a  *"  h 


A15 


%  Inserting  (11)  end  (9)  in  (lO)  we  obtain 


lta[2-(-l)^sin  a]-*4kr, 
ic  jL 


2*3t(aarj)' 


W 


(12) 

•ec§(0j-f)  sec 


_1 

We  nuts  tbei.  u^(P)  is  of  order  k  ‘‘end  u^(P)  16  of  order  k*1. 

Clear  lj  Uj(P)  *111  be  of  order  k’^2.  Bore  u*(P)  is  the  leading  tern 
of  the  field  corresponding  to  the  J-tuply  diffracted  rays.  It  too 
consists  of  e  m  of  tvo  wares.  31nce  Is  of  order  k*^2  it  Is  of 
tke  sane  order  M  the  ...and  tens  in  seen  of  tbs  singly-diffracted  waves. 

We  base  not  computed  these  terns  because  so  far  ee  ere  viable  to  depute  t>» 
anplltulc  eocfflcier.»r.  If  *e  denote  by  u#  the  geoeetrlcal  optics 
field  (l.e.,tbe  Incident  end  reflected  fields)  «e  eey  write  the  solution 
of  the  problem  of  ilffroctloo  by  an  infinite  silt,  with  boundary  condition 


u  •  0  la  the  torn 


u 


(13) 


Although  the  leading  terns  of  the  reenlniag  — Itlply-dl  ffmaetad  waves 
ere  ee  larger  than  teres  onltted  in  (13)  it  1*  interesting  to  note  that 
they  ere  easily  confuted.  In  fact  the  resulting  series  is  a  gas-trie 
eerl*«th—ee  la  easily  eua a»d  (•—  [2<^J ). 


v*  tar*  not  wgww  SnuMy^M  rrmeten  wave  for  the  silt  preblea 


with  the 


condition 


0,  t 


loner  alga  In  (3). 


to  t— 
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Zs  this  c«>e  (d)  vanishes  vher  a  *  x/2,  This  is  to  be  expected,  for  if 
*  plane  uere  travels  toward  a  half-plane  In  a  direction  parallel  to 
tbe  plane,  tbe  incident  plane  vave  itself  satlrfles  the  boundary  condition 
condition  «  C  and  no  diffracted  wave  Is  produced.  If  an  arbitrary 
tnve  u1  is  incident  in  the  sane  direction^  assume  that  the  diffracted 
wive  Is  proportional  to  ,  the  noraal  derivative  of  the  incident  vave 
et  the  edge.  The  proportionality  factor  is  a  new  diffraction  eoafflclnar 
vhieh  can  be  Obt  lined  by  solving  an  appropriate  canonical  problem.  tbl3 
nev  coefficient  and  its  application  are  Riven  in  M. 

Thus  far  ve  have  considered  only  problems  with  straight  edges,  fbr 
e  curved  diffracting  edge.  let  r  denote  distance  along  a  diffracted  ray 
from  the  edge.  Then  the  leading  tens  of  the  diffracted  vave  le  gives  by 


(l«0 


.*  -  s1*, 


(15) 


and,  fron  (7.9), 

»*<«■)  -*»J(0)  )]"^*  (d>»*[r(l«|^)]*^.  (l£) 

Za  CO)  and  (i&),  sl  and  t*  denote  the  phase  aad  amplitude  of  the  Incident 
vave  at  the  point  of  diffraction.  1.  given  by  (7.1?).  If  the  diffracting 
edge  is  the  edge  of  e  thin  screen  aad  the  boundary  condition  on  the  eer 
le  eeO  or  ^  *0,  thsa  (d)  la  given  by  (]).  If,  in  a  aaigHM^houd  of 
the  point  of  dlfftaetloa,  tbe  boundary  is  locally  wedge  mspe-'c,  ihaa  (d) 

U  given  by  (>). 
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To  illustrate  diffraction  by  curved  edges,  we  consider  the  problem 
1  lkx 

of  a  plane  wave,  u  ■  e  ,  normally  incident  upon  a  plane  screen, 
x  •  0,  containing  a  circular  aperture  of  radius  a.  Tbe  geometry  can 
be  visualised  vith  the  aid  of  the  second  figure  of  this  section.  Then 
a  •  0,  and  two  singly  ‘diffracted  rays  pass  through  every  point  Y. 

They  cane  from  the  nearest  and  farthest  points  on  the  edge.  The  angle 
of  incidence  0  ia  everywhere  «/2  and  the  radius  of  curvature  0  of  the 
edge  is  a.  For  both  diffracted  rays  the  angle  6  betveen  the  ray  and 
the  norml  to  the  edge  (which  lies  in  the  plane  of  the  eparture  ) 
satisfies  &  ■  0  -  x/2.  dance  (7.1?)  becomes  t.  ■  -a /sin  3,  Than  (3), 
(7.0)  and  (l’i. 3)  yield  the  "singly  diffracted  field", 

(17) 

UeTj*  iw/U  -J 

“773  I  «,  i  '"I  “7  ■*»  °j3 

Bare  ve  have  added  the  contribution  corresponding  to  the  two  singly* 
diffracted  rays  passing  through  P,  On  the  x*exls,  r^einO^  •  a*  hence 
the  last  Actor  in  (17)  is  infinite.  This  occurs  beoause  the  axis  le 
a  caustic  of  the  diffracted  vavtt.  Since  the  exact  solution  of  the 
problta  is  everywhere  finite,  a  better  asyrgrtotic  eatNuduu  is  required 
in  ths  neighborhood  of  the  axis,  Sue  1  expansion*  art  discussed  in  M. 
If  the  apsrture  ,  instead  of  being  clnr^ar,  is  foroed  by  a  nooth  convex 
curve,  (17)  is  essentially  unchanged,  Ajpin  two  singly  diffracted  ray* 
pass  through  each  point  P,  snaastlng  free  the  nsarsst  and  Aribest  points 


on  the  edge  of  the  aperture .  The  singly-diffracted  field  will  be 
given  by  (17)  If  ve  Interpret  the  angles  and  distances  In  the  obvious 
way.  In  each  tern,  a  oust  be  replaced  by  a  ,  the  radius  of  curvature 
of  the  edge  at  the  point  of  diffraction. 

If  a  place  «*»•  is  normally  incident  upon  a  plane  screen  containing 
an  apvrture  ,  the  edge  of  vhlch  is  on  arbitrary  regular  curve,  the 
diffracted  rays  emanating  from  each  point  of  the  edge  lie  in  a  plane 
perpendicular  to  the  edge.  The  envelope  of  these  planes  Is  s  cylinder 
with  generators  norsml  to  the  plane  of  the  screen.  This  cylinder  is, 
of  course,  a  caustic  surface  of  the  singly. diffracted  wave.  (The  other 
caustic  Is  tbs  sdgs  Ituelf . ).  The  c roes- section  of  the  cylinder  Toned 
by  its  Intersection  with  the  plans  of  the  screen  is  a  curve.  This  curve 
is  the  envelope  of  the  noreli  to  the  edge,  l.e.,the  evclute  of  the  edge. 
Thus  on  every  plane  parallel  to  the  screen  the  caustic  lnteraeuts  the 
plans  in  the  rn&ute  curve,  sad  on  would  expect  to  find  oorreepondlm 
bright  lines  In  tbn  diffraction  patterns  Toned  on  such  planes.  Thee# 
bright  lines  have  been  observed  and  constitute  an  interesting  experimental 
confirmation  of  cur  theory.  When  the  evdUite  lice  within  the  apertore 
curve ,the  lines  ere  naked  by  the  presence  of  the  incident  wave  passim 
through  the  aperture  .  In  such  case*  they  an  ts'srs  easily  observed  vher. 
t;w  aperture  rr.-lscrl  by'tHe  •?'»?11««*rt'*ry  vrwn,  e.  T.^when  the 
circular  apmrture  Is  replaced  ty  a  elreular  disk* 


Al6.  Expansions  containing  exponential  decay  factors  and  fractional  powers 
The  asymptotic  solutions  of  problena  for  the  reduced  vave  equation 
which  ve  have  considered  so  far  have  been  based  on  an  expansion  of  the 
form  (9.1).  However, noro  general  types  of  expansions  have  been  discovered 

wjf  •ijrujtlOlj.C&xiy  O^Cw  SOxuwlO&S  CT  **1C  rCCLUCCm  -~YC 

^  u  +  k2u  -  0,  (1) 


for  a  bosugeneous  mdlua.  In  [7]  >  Frlsdlander  and  Keller  have  mde  a 
syeteeatlc  study  of  asynptotic  solutions  of  (1)  of  the  fox* 


u  ~txp(lks(X).kap(X)) 


(2) 


Here  a  and  1^  are  real  nusbers  and  X^4l>  Although  foraal  solutions 
of  the  type  (2)  exist  for  all  values  of  -  -Mily  the  values  o  •  0  sad  a  ■  | 
have  occurred  fc  actual  problem.  Since  the  eese  a  -  0  reduces  to  the 
expansion  (5.1),  ve  My  restrict  our  attention  hare  to  the  ease  a  ■  j. 

Thus  ve  consider  asymptotic  solutions  of  the  redueod  vave  equation  (1.9) 
of  the  fox*  (1.8)  where 

.  .  (3) 

Since  ve  have  ahovr  that  a  satisfies  (1,9)  ve  my  Insert  (3)  in  (1.?)  to 
<Wa»1n 


00 


-ka[C*)2-nV21k^3  v*.^  ♦  Ik^f^vwAs] 

♦  k2^3  (?*)W^3(2W.*ne^)«*  -  0. 
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Proa  the  fora  of  (b)  it  It  clear  that  we  nny  expect  that  v  will 
have  an  expansion  In  reciprocal  powers  of  klyl^.  Thus  we  sat 


•r  S  m  fvYV  ^  «  X  _  9 

(1) 

/  a  Z_  “  '  ' 

n»0  n 

where  s  a  0  fur  a  •  -1,-2,.,,  .  If  vs  insert  (5)  In 

(’4)  and  collect 

powers  of  y,  obtain  the  equations 

C*)2  «  n2, 

(6) 

a  0, 

(?) 

and 

^a’^  4  *  V 

(8) 

rm  ‘  ’  ^^n./^4*®.^  41  * 

Ve  note  that  (6)  It  tha  familiar  olconal  equation.  It  follows  that 
t*<a  rain  ftaturet  of  our  earlier  expansion,  lae.,the  rays  and  save- fronts, 
are  pcttsmd  in  the  new  expansion.  (7)  aerely  asserts  that  the  surfaces 
p  m  const  ere  crtboconal  to  the  wav*. fronts  t  •  const.,  !*•* * 

p  a  const,  on  each  ray.  (10) 

Pgr  n  a  0,  r^aO  and  (8)  Is  Identical  to  tbs  sero-ordsr  transport  tqtstlon. 
for  arbitrary  a  (8)  oaa  be  written  In  the  fora 
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dz 


2n  4  *  ^  *  r  . 
do  m  01 


(n) 


Sere  a  denotes  arc  length  along  a  ray.  By  conpariaon  with  (3.14)  we 

easily  obtain  the  solution  c?  the  ordinary  differential  equation  (ll), 

(12) 

t(«’)n(e')lV2  rja.) 


,  .  ,  .  rg(roH°o:  r“  1  f  t(o’)n(a')l 

,(*)-^(ff0)  [pT-iirai — J  45  4o  [TTTO1  j 


do* 


for  a  hoaogenaoua  aediua,  n  -  const.,  and  (12)  bee ones  (See  section  A4) 


,  [ fi  r<v*,)(po«,nt/2 

*.(9) -l*(%)  [TS^fe^j  *s  l  j  r. 


(13) 

(a'Jde*. 


Our  new  expansion,  which  takes  the  funs 


u'Vexp  Uks(x).k^"V*))  T**  *_(X)h  J  , 


(14) 


will  also  be  oalled  a  It  will  be  required  shortly  in  our  discussion 

of  liffraetion  by  aeooth  bodies. 


i 

i 

i 

j 


A17.  The  surface  eiconal  equation  and  rurthca  raw. 

In  prepare t lor  **r  o»r  study  of  diffraction  by  snooth  bodies  we 
consider  now  the  Initial  value  problaa  for  the  eiconal  equation  on  a 


surfhce.  Ve  are  concerned  with  a  function  s  defined  only  on  »  »urf«ce  3, 

end  with  initial  values  prescribed  w*  a  curve  which  lies  cn  that  surface. 

Let  X  ■  x(t. ,  t  )  be  a  pa react ric  equation  for  the  regular  surface,  S. 
x  2 

following  the  custceary  notation  of  the  differential  geanetry  of  surfaces, 
we  introduce  the  surface  target &  vectors 


v  *  ,  «  . 

X1  ■  ^  '  X2“^1 


»rvt  the  metric  coefficients 


-  VX«*  *’4  •  **  W 

Vs  also  introduce  the  lnverae  (g1J)  of  th#  setrlx  (gjj).  Th*R» 


course 


*1 3  " 


Xa  (j)  sad  subsequent  equations  we  asploy  the  sua»tlon  eorvention  for 
reposted  Indices  over  the  values  i,2.  Is  the  Krosnscher  cy*ol* 

df 

P Br  any  function  ftr^Vj)  of  tho  ourfhoe  pamsMtcra,  Ut  fj  •  ^  . 
Tbs  surface  gradient  of  f  la  defined  by 

*  .  /l  rt  xk.  <*> 


ts*  a—  that 


(^)  agrees  with  ths  us«l  deflaltloae  of  tho  gradlss*  ws 


set  dX  -  X¥dvv  sad  ctessrws  that 


J 

* 

I 


A’ 


(?> 

Sr.u .  .klf,WT.  •  ■  'iV’v  ■  V’,  ’  "• 


X-  r.  « i -  »k«  Inin  nf 

- ...* -  .a.  *»  -j.  j-  — . . . 

refraction  nfr^  fg)  •  Tg)J,  the  ourfacn  elcooal  equation  cen  be 

written  lh  the  equivalent  fores, 

„  (6) 

C^)8  ■  n2|  h^*l,fc8,Tl  t2^  B  *i^TllT2^*i,j*n*^Ti,T2^  " 


la  order  to  eolve  the  flret  order  partial  differential  equation 
(6),  we  Introduce  the  cheneterlstle  curvet  [^(e),  ra(o)]  which  ere 
determined  by  the  eolutlone  of  the  cbereoterletlc  equation*  (taelltoa'e 
etuetlcee) 


*1  •  3  3^  -k«1JV 

•l  ■  *3  it  *  ’  3  C<A  Vj*fa2)J- 


(7) 

(*> 


the  dot  deaotee  differ eutUUot.  with  respect  to  the 
(6), (7)#  (*)  U*ly  tint 


.a 

x 


*. 


(9) 


*  "  *thTl’h  *  ^  ^  xfjf***  •  ^8IkV*^< 
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Banco  *e  vey  identify  the  parameter  o  vita  *rclmgi.h  «luug  the  surface 
cvrv»»  X  -  X(o)  -  X|j1(a)>  Tg(o)]  by  eettlng 

*  «  £  .  (10) 

These  curve#  vlU  be  railed  aurface  rtw.  Wa  not*  that  the  equation 

*  •  *i  *i  ■  “1J*A  • :  *  <u» 

lspliaa  that  the  surface  ray*  are  everywhere  orthoemel  to  the  aurface 
♦ave- fronts  aCt^tg)  •  coast,  fraa  (7),  (6)  and  (10) 

a  •  s^vi  ■  IfM  a^Sj  •  la*  •  n,  (u) 

benev 

•[X(c)J  -  *Dt(c0)]  ♦ 

(u)  provides  Uat  solution  of  the  surface  eicoeal  equetlnn  (6),  ones 
initial  valuta  are  epee  1  fled  , 

Wa  aaaoaa  that  the  Initial  values  are  fives  on  a  curve  oa  the 
surface,  (t^Vg)  -  (V^Ol),  VjCttf],  vhere  T«  la  as  arelanftb  foiasatar. 
Haw  the  Initial  data  take  tbs  fora 

spCWl  •  e°fo).  (U) 

lore  «°CU  it  a  fives  function.  OUferesUwu  of  (IV)  yields 


/ 


»[*(••  )J**\ 


(13) 


(13) 


.  1  da° 

001  p-s  w 


UaM  A  4.  »Ka  mwm  ■  ««•*».«•  *•?  «~4  f)w>  initial  curve.  If 

«m  AAituiM  that  «1  <  ■■  <  1  tbu&  (lG)  implies  that  rt  every  point 

on  the  initial  curve,  cm  rurface  ray  i>  outgo!  **  fr<»  that  curve  on 
each  aide  of  it,  'Awm  surface  raya,  together  with  (13 )  provide  the 
outgoing  aolutien  of  the  initial  value  jroblen  for  the  surface  el  canal 
equation, 

wv  rwoall  tint  the  r-ye  associated.  with  the  elconal  equation 
(l.ll*)  becoae  etraieJit  linoa  in  the  caae  n  •  const,  Iheae  straight  liaee 
are,  of  course,  aoodeeics,  or  ahorteet  paths  between  two  points  in  apace. 
We  will  now  prove  that  for  the  ease  n  ■  nonet,  the  surface  rape  defined 
by  (7 )  and  (8)  are  geodesics  of  the  surface  5,  The  proof  will  occupy 
the  renainder  cf  tfciw  section. 

If  n  •  conet.  (7)  and  (8)  take  the  row 

\  •  i  «iJ,j  *  4j  ’  *  fe  (^v,j  V/  {l7) 

We  replace  thic  eyeten  of  four  fir*t*order  ordinary  dlffneotlal 
ipatloae  by  a  eyeten  of  two  eenani-oicer  eeuatioM,  We  first  note  that 
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•  •  1  »ik-  «  „<5y  • 

Vj  "’2  s  *k6  av  * 


Next  ve  introduce  the  Christ off  el  symbol  rr)  defined  by 
f^r1!  ■  i  an<(sr  ).  *  (a.  ).  -  («..)  ). 


end  (19)  Md  (20)  field 


t^1)  Vj  "  T2  f’ik8'3wen,^*ja,i  +  ^iJj  ’  Vv  • 


This  novation  can  be  simplified  by  usin'  *n  identity  obtained  by  differentiating 


It'h  *ij  ’  •«kl(«iJ)v 


Tbe  result  le 


(ViVj  •  •  £  (.“."va  •  «*vvA  • 

.  -  -S  («“(«"),  •  •  «"A )  vv  <«i 

ft) 

rmn  (ifl)  end  (23)  ee  easily  obtain  the  second- order  eystea  of  differential 


eguatiom 


v;  ♦  (V)  ^  -  o. 


T*T» a3.  J.1* » •; OS V. 


*  See»  e,in  Stoke?,  W»i 
notes,  (1955),  pegs  T-5,  elation  V.b. 


[  A  i  vW 


Ifcw  fork  University  Lecture 


f  J 

* 

(0 
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A  18.  Diffraction  by  smooth  objects* 

In  this  aectlon  w«  »h*ll  derive  a  gancial  formula  for  the  diffracted  vurt 
which  is  produced  when  a  vave  u1  le  incident  on  a  smooth  surface  8  In  such  a 
way  that  some  of  the  incident  rays  era  tangent  to  8  along  a  curve  C.  In  this  cast 
there  Is  a  shadow  region  which  is  not  penetrated  by  any  of  the  ordinary  rays  of 
geo&etrical  optica.  The  shadow  region  Is  separated  from  the  region  reached  by 
incident  and  reflected  rays  by  n  surface  called  the  shadow  boundary.  The  tangent 
rays,  beyond  their  points  of  tangeney,  lie  on  the  shadow  boundary.  For  slrruiclty. 
we  shall  assume  that  8  la  a  boundary  rather  than  an  Interface.  Thus  we  shall 
avoid  the  additional  complications  of  transmitted  waves.  As  In  all  of  our  con¬ 
sideration*  ,  the  following  construction  will  involve  certain  apparently  erlbtrary 
prescriptions.  These  prescriptions  wars  discovered  by  examining  the  asymptotic 
expansion  of  exact  solutions  of  boundary  value  problems  for  the  reduced  *«v« 
equation.  They  will  be  further  trifled  by  the  boundary  layer  theory.  However , 
general  proofs  of  the  validity  of  the  formula  have  not  yet  been  given. 

In  order  to  derive  the  formula  for  the  diffracted  wave  u5  we  firet  con¬ 
struct  e  surface  vave  (or  creeping  wave)  ue  vhten  Is  defined  only  on  the  surface 
8.  The  curve  C  eets  as  the  (secondary)  source  of  the  surface  wave,  which  Is  e*. 
cited  by  the  Incident  w*vt>  u  .  ue  Is  defined  only  uu  the  "dart"  side  rf  C,  l.e., 
on  the  portico  of  8  adjacent  to  the  shadow  region.  On  this  portion  of  9,  the 
phene  se  of  the  surfeew  wave  satvsrits  the  surface  eleonal  equation  (17.6)  with 
Initial  conditions  given  by 

s  -  s  on  C. 

of  the  eateriel  In  this  sect  too  Is  adapted  from  (36). 


(1) 
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It  follows  cosily  from  (l)  that  at  each  point  0^  on  C  the  surface  ray  emanating 
from  that  point  la  tangent  to  the  Incident  ray  (which  la  tcngcnt  to  S  at  Q^). 

If  P^  la  any  other  point  ou  the  aurface  ray  emanating  from  Q^,  we  see  free  (17.13) 

^  a^)-. (2) 


e. 


Here  the  variable  of  integration  0  la  arclength  along  the  aurface  ray. 

Before  finding  the  amplitude  of  the  aurface  wave  ve  will  begin  the  con¬ 
struction  of  the  diffracted  wave.  The  "dart"  aurface  of  3  acta  aa  the  (secondary) 
source  exciting  ud.  The  phase  sd  of  the  diffracted  wave  eatlwflea  the  elconal 
equation  (l.lb)  with  Initial  data  given  by 

ad  ■  ac  on  8.  (3) 

We  aee  from  section  6  that  ed  la  th«  solution  of  *  "characteristic  Initial  value 
problem"  for  the  elconal  aquation,  and  that  at  every  point  P^  on  8  the  diffracted 
ray  eannetlng  from  7^  la  tengant  not  only  to  8  but  also  to  the  tun  ce  ray  passing 
through  7|>  Portion*  of  the  Incident,  eurfocr.  end  diffracted  rays  ore  sketched 
la  the  following  figure. 


The  rays  may  hr  described  a*  follows:  the  incident  rur  which  la  tangent  to  8  at 
^  splits  Into  two  brvnrhea.  o*>  branch  (not  shown  tn  the  figure)  continues  along 
the  shadow  tabluy)  the  other  branch  Is  the  surface  ray.  At  every  point  ?j  (on¬ 
ly  one  point  la  shown)  on  Its  path  the  aurface  ray  splits  Into  two  branches.  Oik 
braoch(not  shown)  continues  along  the  surface;  t**»  other  branch  1*  the  diffracted 
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ray  emanating  from  p  .  Prom  (6. '2)  the  phase  of  the  diffracted  vave  la  given  by 

a  r ^  *  |'P»  r  ? 

s1(P)  -  aC(P1)  ♦  J  wlo  ♦  J  nd a  *■  J  ndo.  (4) 

J  P1  J  ®1  P1 

Q 

The  leading  term  of  the  surface  vave  is  given  by  uc~  e1  tQ.  In  order  to 
construct  the  ooplltudn  *£,  ve  consider  the  vldth  dv{o)  of  an  inf lr.it ec Ira  1  strip 
cv  mi.  face  rays  st  the  point  o  on  a  civen  surface  ray.  The  'fenergy  flux" through 
such  a  strip  lc  proportional  to  n(<?)  [*£(<*)  ^dv^).  At  th?  point  o*d"  v*>  nssui.ie 
that  the  n>iv  <-  W  due  to  energy  lost  to  the  diffracted  rays  uhleh  emanate 

free  ih#  surfsee  reys  in  the  lnvervai  do,  and  thax  the  energy  lorn  is  propomunax 

c  2 

to  n(<0)  and  to  the  area  slsnent  tivoo .  Tnus 

atn^)2  dv'  -aa  n(t®)2  dvdo.  (5) 

o  o 

The  decay  exponent  cKc)  depend*  on  local  properties  of  the  surface,  the  medium,  and 
U'.e  field.  Integration  of  (5)  yields 

<W  *  «■»  (-]“  «!••>*.•).  CM 

lie  aw.  oseune  thst  the  amplitude  cf  *he  surfsee  veve  et  Q,  Is  proportional 
to  the  amplitude  of  the  incident  vave  at  thot  point, 

4^-  *  (7) 

■ere  d(^)  V#  «  4  f fraction  coefficient.  Horn  (f.iy)  «*  onteln  the  formula  for 
the  aapiifcuoe  **{p)  of  the  diffracted  >avs  at  a  distance  3  alocg  the  diffracted 
toy  fro*  the  po-'t 


g,  ,  ~d,  ,r  .  o(on? 

i0(s)  -  si#  7  -3351  sjyjj 


The  guantitice  0^,  v,  d*^.  ana  are  defined  la  section  T.  esawse  that  the 


amplitude  or  ‘  lie  diffracted  rfayo  i;,  proportional  to  that  of  the  surface  wave  at 


the  point  , * 


za( 

0 


-  ;f-'-;(r1!zo(P1)  .  (9) 

The  diffraction  coefficients  d(Q^i  and  d(P^)  are  assumed  to  he  the  same  function 
of  the  properties  of  the  surface,  medium,  and  the  field  at  the  respective  points 
4.  and  r. .  This  assumption  is  h»«~t  on  the  reciprocity  principle*-*  source  at 


Q  produces  the  same  field  at  P  that  a  source  at  P  produces  at  Q.  We  shall  see  that 

the  values  of  the  diffraction  coefficients  and  the  decay  exponent  depend  on  ..ie 

boundary  condition. 

Fran  (6  -  9)  we  now  obtain 

*  r }  \ 

i  \  J 


z^Pl-k'-dlPj  »d(H,  )r*(Q, ) 


)  nfi)  liO  d#, 

o^r  c2 #in '  sttf* 


Then  ihv  l«alin£  Lera  of  the  diffracted  field  is  Riven  by  (I*),  (10),  and 


J(P) 


iU*tP) 

* 


;n) 


TV**'*  '-muttons  were  derived  uni  nr  the  surface  wave  uc  -  e**#  *C,  which  no  lower 

W 

appears  in  Uw  i-eeult.  Tie*  iuIiw  wav*  is  not  an  ssyuptotic  cup; .tier  ?'  it* 
true  solution  st  the  boundary.  (This  is  especially  clear  if  the  boundary  condition 
lit'  1|  It  is  nerely  «  convenient  Intermediate  step  tc-  the  •'•scriptiau  of  viw 
Ulffractcc  vcvc.  The  latter  l*  sinrular  e»  the  b^x-m'iae  the  diffracted  revs 


have  a  caw. ic  ''.ere. 


TUr  cneti  rtjor.  sc  ffcr  is  not  mill*  complete.  Act>*lly  the  diffracted  wave 

d 

conal.-te  or  a  n«t«r  of  »o>a  u^{  J  *  l,?,  ...  of  which  we  have  construct*-'.  ‘dUir  ore 
In  (11).  hich  tode  nae  Its  own  diffraction  coefficient  4.  and  4~c»y  expos.- nt  »t. 

Thus  th>*  ilffrac^ed  f'.vl-l  ll  by 

- - T - 

The  factor  a  '  lr.  (9)  Is  required  in  order  that  iiP(;  slould  be  dinensloniese. 
iSe*  It  p<*  -ion 


I  *1 

^(4,)  exp  ^  ik  [*3^(0^)  +  J  ndo  +  j 

L  4,  ''P 


au(u  )  n(Q.) 
dw;p1)  "n(p)  p?  sir‘  7 


dG,dS. 


LdS.-j^  ■r-i  i 

fefj  z k  “WW 


-  J  a.(o)dc 
%  '  - 


for  the  case  oi  a  homogeneous  medium,  n  ib  eoustiiYt  and  (7*21)  and  (12)  j-!?1 
ud(P)  ~  z*(Qj  exp  ^ ik  12(4.  )  +  nr  +  no  1  \>  ^7T~!~^  1 

°  L  L  JJ  l“  '*!'  "*  j 


x  k*^  dj(P1  )d1  ( *a  )  exp  -  j'  <i. . (o )ilu  j- 

.1  '  L  \  '  J 


Here  a  Is  the  distance  from  P,  to  P  and  t  is  the  distance  fran  ^  to  P  along 
the  surface  geodesic.  Values  of  d^  and  will  he  given  in  (19.13)  and  (19.1'i) 


of  Gection  19. 


■e  use  (7.21:)  instead  of  (7.19)  deriving  (fl)  then  (12)  takeo  the 


ud(P)  ~  /((^  exp  aids  (^)  +  /  ndo  ♦  /  nd-.j  S 

L  \  V,  J 


iw(Q. )  n(d.)  da(P,  )1* 


.A 


2 jk  exM  "j 


r1  .  \ 

i  a  ( o  /da  >  . 


We  see  from  (7.23 )  t.'iat 


.  am  _ q&H-L  (15) 

da(P)  P’-P1  coda(P)  • 

Here  f  is  a  point  on  the  diffracted  ray  Joining  and  P  and  cQ  denotes  the 
distance  from  to  P1  along  this  vay. 


-75- 


A1  ’J 


A  IS-  diffraction  tv  a  circular  cyi i nder 

From  (18.7)  urd  (18.9)  we  see  that  the  diffraction  coefficient  <i^ 
is  dimensionless  and  (lO. 5)  shows  that  the  decay  exponent  has  the  dimension  of 
a  reciprocal  length.  The  diffraction  coefficient  must  depend  on  k  because 
we  expect  it  to  vanish  for  k  -*•».  Thus  d^  must  be  a  function  of  he  where  a 
is;  a  Ikjjos,;;.  jssiir-:  tv.-.*  far  •  hcmc.-enec1"®  "wiium  els  the  radius  of 


curvature  of  the  normal  section  of  0  in  the  ray  direction.  Wc  aloe  assume  that 
ttj  depends  only  on  k  and  a.  Then  dj  and  can  be  obtained  from  the  asymptctl . 
expansion  of  the  exact  solution  of  a  problem  with  some  simple  surface  8.  In  this 
section  we  shall  find  the  expression  for  the  field  produced  by  a  line  source 
which  is  parallel  to  a  circular  cylinder  of  radius  c  in  a  medium  vith  index 
of  refraction  n  ■  1.  Comparison  with  the  exact  solution  will  yield  the 
coefficients  dj  and  a.j.  The  problem, which  1*  two-dimensional, is  Illustrated 


Let  r  denote  distance  from  the  source  point  Q,which  is  located 
at.  the  point  with  polar  co-ordinates  (  P,  0).  We  take  the  Incident  wave, 
produced  by  this  source  to  ba  (compare  section  9) 


fo(1V) 


(1) 


*  The  meteriel  in  thie  section  is  besed  on  (S9) . 


-76- 


h 


The  surface  rays,  which  are  geodesics  on  the  cylinder  are  clearly  arcs  of  the 
circles  which  generate  the  cylinder,  and  in  (1C. 13)  it  is  clear  that 

dv(<L) 

dw(p^)  *  1»  4111  ^2  *  “*  The  »®sunptions  made  above  imply  that  a^cr)  ■  const, 
and  d.(P. )  u  d.(Q. ).  Furthermore,  since  Q,  is  a  point  of  tangency  or  a  ray 
firca  Q,  we  see  that  at  ,  r  “  (p2- a2)1^2.  Similarly  «  ■  (r2-*2)^2. 

Thus  (18.13)  yields 


ud(r,$ 


*) -[sskV-2)1/8^-.8)!/2)-!/2  ^[(pV)^)1^} 


*  y~  dj2  exp  {(ik-aj)  T  ).  (2) 

3 

I 

3 

Equation  (2)  gives  the  field  on  a  ray  from  Q  to  P  having  an  are  of 
length  t  on  the  cylinder.  For  the  ray  Q  Q^P^P,  t  ■  to  where 


t0  ■  a®  -  a  coe*1(a/a)  -  a  co#*1(a/r). 


13) 


»it  all  rays  which  are  tangent  at  encircle  the  cylinder  n  times,  end  leave 
at  Fj,  also  contribute  to  the  diffracted  field.  For  these  rays,  t  -  Tn  where 


t  »  t  +  2n«a. 
n  o 


(h) 


We  note  that 

n«0  *■' 


(3) 
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Therefore  we  may  insert  {!»)  in  (2)  and  sun  over  n.  This  yields  the  field 
contribution 

!6) 

*  y~  dj2  expl(ih-ttj)  Tj  [l-expCSitaflk-a^))]’1  . 

At  every  point  P  there  is  also  a  contribution  ud  corresponding 

tc  rays  which  encircle  the  cylinder  n  times  In  the  apposite  direction  end 

leave  at  P  .  can  be  obtained  by  replacing  6  by  2s-®  in  (6).  Then  the 

total  diffracted  field  ud  »  u?  +  ud  la  given  by 

^  2  (7) 

ud(r,®)*~  (6s)*-k”\r8-a2)*1^  (p2-*2)*1^  exp{ik[(e2-a2)1^2+(r2-a8)1/raJ+l*/U) 

x)~*  dj2  [l-exp(2«a(ilt-ej)fl*1[exp{(ik-aj)a®)+exp{(ik-Oj)a(2s-®)}] 

3 

x  exp  (-(ik*aj)a[coa'1(o/o)+cos*1(a/r)3 ). 

Except  for  the  coefficients  d^  and  a^,(7)  is  an  explicit  formula 
for  the  leading  term  of  the  diffracted  field.  In  the  uiiulow  region,  this  la 
the  only  field.  In  the  "lit  region"  it  must  be  added  to  the  incident  and 
reflected  fields.  The  coefficients  d,.f  and  depend,  of  course,  on  the 
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A19 


boundary  condition  specified  on  the  cylinder  r  ■  a.  We  will  take  thia  J 

condition  to  be  the  impedance  boundary  condition  *  ii«tt  ■  0.  Hare  >i>>  ; 

« 

constant*  ' 

r  ,  3 

In  ]29J  the  above  problem  is  solved  exactly  by  separation  of  « 

$ 

variables  and  expanded  asymptotically  for  large  ka.  The  result  agrees  exactly  „ 

with  (7)  if  we  set 


and 


■  ’'W6  ^  V 

dj  ■  .*1/s  c  {*’(*))}  *♦ 


(<J5 

(9) 


*th 


Here  q^  io  the  J  solution  of  the  equation 


r$  -  •  *1/s<V/3«. 


A  (qJ 

and  A(x)  is  the  Airy  function 


(10) 


m 

AM  .j  coa(r^-xr)dT. 


(ID 


A  (x)  denotes  dA/dx. 

H»d  we  chosen  a  constant  index  of  refraction  a,  other  than  n  •  1  it 
is  dear  from  the  form  of  the  reduced  w  ve  equation  and  the  impedance  boundary 


« 
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condition  that  (8-10)  would  ha  modified  by  replacing  k  by  kn  and  s  by  z/n. 

In  order  to  determine  the  diffraction  coefficient  d^Jx)  and  the 

decay  exponent  a  (X)  (at  a  point  X  on  the  boundary  aurface)  for  the  case  of  an 
J 

inhomogeneous  median,  comparison  with  other  exact  solutions  (See  [3^J  and  ) 
indicates  that  we  must  replace  l/a  in  (8-10)  by  the  "relative  curvaturo" 

[l/a(X)  +  *(X)]of  the  surface  and  the  diffracted  ray  emanating  from  the  point  X  . 
Here  e(X)  is  the  radiua  of  curvature  of  the  normal  section  of  8  at  the  point 
X  in  the  direction  of  the  surface  ray  at  that  point,  and  x(x)  la  the  curvature 
of  the  diffracted  ray  at  X.  Wa  now  make  the  replacementa. 

i  *  *(X)  ,  k  -*lm(X),  >  -**/n(X)  (12) 

in  (0-10).  The  result,  aftar  soma  simplification  of  (9),  it 

dj(X)  -  a'1  Aj  [6’1  ta(X)]l/j  [a*X(X)  ♦  *(X)]8/3,  (13) 

d,(X)  -  -U/Zk  6'1/6  Zml\^tka(X)f/6  (a*1(X)4<c(X)]‘l/6  (Ik) 
J 

Bare  q.  ia  tha  Jth  solution  of  the  equation 
J 

„  .5*1/6  .k1^  |6n2(X)[a'1(X)4<(X)]  J  *l/3  .  (13) 

i 


v 

w 


V 
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The  boundary  condition  used  in  det* **ndnin|  (13*19  ie  the  impedance  boundary  condition 
^  ♦  ik*u  »  0  ,  on  S. 

Bare  ^  denotes  the  normal  derivative.  If  s  ia  net  a  constant  on  8  it  must  be 
dv 

replaced  by  s(X)  lvi  (1>). 

2.  ^ 

Vj  tee  frcn  (13)  that  a,  is  of  order  k  1  .  Therefore  the  formulas 
(18 .12-1*1 )  for  the  leading  term  of  ud  agree  vith  the  general  fora  of  expunaioa 
studied  in  section  15.  Lower  order  terms  in  the  expansion  for  u  can,  in 
principle,  be  obtained  by  boundary  layer  methods 


ttany  isterestim  iwatures  of  the  foregoing  theory  can  be  illustrated 


by  considering  problems  in  which  the  index  of  refraction  is  constant  on  planes, l.e., 
a  function  of  a  tingle  cartesian  co-ordinate  x,  for  in  this  case  the  ray  equations 
can  be  integrated  explicitly.  Vb  consider  s  problem  with  a  plane  boundary  at 
x  ■  xb  and  with  an  index  of  refraction  n(x)  which  increases  aonotonioally  far 

X^<  X  ,  At  a  ■  we  Impose  the  impedance  boundary  condition,  vtth  constant  z, 


♦  iksu  ■  0  .  M 

A  line  source,  iwrpendimilar  to  the  plane  of  the  fpllawing  figure,  intersects 
thlaplen*  at  the  point  (xQ,  0). 


*Thia  is  strictly  trus  only  in  the  esses  s  ■  0  end  s  •  •.  Otherwise  we  see  rrom  (l?)  *.he«- 
1  ■  a  function  of  k  aid  the  k  dependence  of  otj  is  obscured.  If,  however,  we  set 
z  .  k'^^p.  where  *o  Is  indepe»ew:i'  of  k  then  is  independent  of  k  and  the  above  state- 
ment  Is  trj«  for  all  iQ. 

**The  materiel  tr.  tv-in  ,»ct.ion  le  based  or  f3*1* 
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,?  (*>/) 


Q^vd 


!  ?o\»X 

s(Ci 

t  / 


(xI^T  (*<>.°) 

Aa  la  Section  9  we  shall  characterise  the  source  by  the  inhoaogeneoua  reduced 
wave  equation  +  kV(x)u  -  -C(x-x0)6{y).  It  suffices  to  confine  our  attention 
to  the  const! « :tlon  of  the  *’sld  In  the  upper  half-plane  y  >  0. 

If  we  set  1  ■  1  and  danote  the  parameter  la  (2- **»5 )  tor  -»  then 
the  ray  equations  take  the  fora 


$  ■  Vt  >•  $  -  «.  *  <8>*  ’  »*•  <al 

at 

2  2 

to  Integrate  (2)  It  Is  convenient  to  set  v  ■  J*.  Then  ^(f-)  • 

nonce  via2-  a8,  where  a  la  aa  arbitrary  constant,  the  last  two  aquations  la  (2) 

now  laply  that  ■  ^a  .  Since  ^  •  •(>*•  It  fallwa  that 


d* 


(rfTPJ** 


,(.)  •  r(»0l  *  /  ■ 


J 

-V 


1 
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The  v*y*  einanatinc  fra  i  the  eource  vlll  be  called  "incident" 
rays.  Let  tan  a  be  the  slope  of  such  a  rty  at  the  source.  Then  it  follows  from 
(3)  that 


a  »  n(xQ)  sin  a. 


(5) 


If|?/I  <  the  inellient  rays  proceed  to  the  i^^iand  are  given  by 

.  (€) 

1  • • ?  I  ‘  t  • 


Hare  x<  and  x>  denote  respectively  the  eneller  and  larger  of  x  anl  xQ. 

For  ie/».  <  a  <  a  some  Incident  <>iys  hit  tbs  boundary  iml  «.r«  reflected 
’  Mle  the  others  Uecoae  vertical  and  are  turned  b.  eh  to  the  ri&t  before  hitting 
the  boundary.  '*«  eee  frow  (J)  ini  (*.)  that  cuch  te.-oltur  points  occur  "her 
n(x)  »  a  •  n(«0)  sin  a,  Eoyooi  the  tvruinc  point  au~.li  r*yi  -ill  he  r tiled 
refreeted  it  ye.  A  articular  rry,  -nth  C  ■  is  toncent  to  tha  boundrryj  i.e., 
its  toiT.li i  poii.u  is  at  x  -•  Time  la  determined  by  the  squctlcn  1 

-  n(XJ  tin  (7)  j 

This  "Haiti  r/r  ray"  la  UluetrataA  In  the  figure.  Its  continuation  lias  on  the  1 

i 

ahodew  boundary,  Irxldant  r-yn  vltb  </;  <  u<  ^  produce  refrooted  rays  at  their 

< 

turnlnr,  point-.,  Incident  mys  with  a  ■'  s  <  n  produce  rcrltcted  rays  ot  the 

•  •  i 

boundciry,  I  >  e'.dltlon  the  llnltlug  rvy  produces  r  wrf.es  ray  on  the  boundary 
eud  dlffrr.eted  .tys  in  tiw  whsdoe  region. 


In  order  to  csleuL.te  the  phuse  «v.  an  Incident  my  we  use  (t.ll), 
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Thus  we  obtain 


t  x>  „ 

r. 2,.,.u«  ..  r 


Here  we  have  used  the  identity,  dx  -  t(n2  -  a2)*dt,  which  was  derived  above 


equation  (?).  (6)  snd  (8)  now  yield 


pi  “  J 


>..«)*  4../  7^3  ■•>•*/(»'  •=)* 
(n  -  a  )  - 


dx  .  (9) 


The  wave  produced  by  an  isotropic  line  source  was  derived  in  Section  9- 
Therefore  fro*  (9.10)  we  nay  conclude  that 


*o  "  e 


^  [ran  s&r]- 


order  to  calculate  ve  see  fra#  (5)  w*  (8) 


±L  .  iiX  £  .  fn2(*  )  -  a2]'  f  j)  ^ .  (u> 

da  to  to  L  0  J  w/  /nz  aZ)?  * 


and  fron  the  figure  we  see  that 


J  . 

eos  ?  ■  (1  ♦  tan  > ) 


i/!  •  i-tfr-* 


(n?-a2)3.  (12) 


Ain'**  dw  •  dv  rno  ;  It  follows  that 


8  -  ngcosy  -  [n2-  [n^-s2]  J 


By  Insertin'.  (13)  j.n  (10)  ue  uttein 


4  ■  -w‘  [»X>  -  ‘t*  f?’-  ‘t1'*  hfev’l 

Thus  the  leading  ter*  of  the  Incident  -.savt  Is  given  by  (d)>  (it)  end 


:  iks1  i 

u  ~  e  I0  . 


For  <  o  <  *  the  Incident  roy  hits  the  bound,  ry.  The  corresponding 
reflected  rey  Is  obtained  by  reflecting  the  Incident  rey  acmes  the  borlsonfeal 
Una  y  -  ^(x^),  Theref^'  •  It  la  'lven  by 


-  Jf  > 

>  t  J 


.  <l£i 


y  •  /(«)  -  iyt(\)  -yt(x)  -  *j^  J  * 

The  pbose  on  the  reflected  rey  la  ootelnad  by  an  argesant  stellar  to  that 
which  led  to  (9),  The  result  la 

x  (IT) 

/  •  •*(*>)  ♦  ejVVOs,)]  ♦/  (»*-  -2)172  to 

•  •'*[(*(}  <»■’-  •*»*  “ 


In  order  to  detendne  the  reflect**  « solitude,  m  first  ««c  (10.6). 


Thia  yields 


«  "(O  cot  0  -  a 

•o(S>  *  "oK**  r  *  nfit^)  cos  & T7 
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Here  0  is  the  angle  of  incidence  (or  angle  of  reflection).  Aa  in  (12)  we  see 
that 


coa 


0  ■  n'1(xb)  [^(x^J-a^J1^2, 


(19) 


lwnce  the  reflection  coefficient  r  ia  given  by 
r^/-  *  .3,1/2  .  . 

r  ■  i^Z  'w1*  ■  • 

Next  vo  use  ().7)  and  (3.6)  to  obtain 

r  r,  ,  f»(Vdv(S>lV* 

*o  ■  ‘oK*  [  gxTfrftT  J 

and  then  (31),  (lfl),  and  (it)  yield 

i  ■  .•w"  faj-Y*  &V  ^ 

„r»* /*•.!*_  t*|) 

L  \  h1.  ,l)m  L 


(*0) 


(21) 


(22) 


-1/2  fn^)  dw(^) 


1/2 


Us) 


At  la  (13)  we  can  ahow  that 

n(x)  •  [n2.  a^*’ 

hence  (*fc) 

r° f  fr°  A„'«»  T 

sRtesr  ■  prrp. 


ircrrtiii*  ' 24 )  m  we  obtain 


in/U .  2,  .  2. 

re  fn  (x  )-«  ) 


-l/1*  ,  2,  .  2 

'  fn  fx) -■  )  ' 


The  1  re^ir-g  teiT  of  the  reflected  wave  ta  gtvon  by  (17),  (2J-),  and 


(26) 


In  order  to  dcUmrit-e  the  refracted  vavt  vc  eanaldcr  tne  Incident 
rv:  •*tth  start in*;  tnfjc  tn  tl«-  inter-*  >  */-  £  ■'t*  Th*  *t)w  xQ  or 

*  a*  >he  turning  l“‘tn*  J*  1 t tee  hy 

r.  (x^)  •  «  •  n(«w)  cln  *  .  (^i) 

Jfp  u»  the  turn  In#  point,  (6),  (7).  *r*t  (1*)  a  tv  valid,  t*;‘.  (1*;)  in 
inielenal**!*  »i  i  .  Therefor*  v*  wist  deumlr*  the  Unit 


..  I 


XU  UU  tHlJL*  jCl  £,  *  rv- f 


\  r.A  *  t.h*  ojcnancion 


Then 


(29) 


n2(x)  ■  ac{,xQ)  ♦  cnn*(x ^(n-x^..)  ♦  ...  •  o'  ■*  *-'t  ♦  ....  •  Tan’ (x^). 


.  X  -X 

o  * 


X  •  * ...]  |  /  [•^u  3/2+ 

.V*r  ,  i  "l*o*  *cr  90?  n<*o> 

L  J*  *  a'fx,,) 


Si2  Urn 

•*T*  *0 


(50) 


Tw  If  wo  let  x  *  Jtg  in  (it)  tte  raault  la 


(31) 


The  nuontlty  *  tea  an  Interacting  gec—trto  Interpretation.  To  M  thle  *e 
writ*  the  ray  tratlum  (2,6)  In  tte  vutor  for* 

it  •  • 


•’  K  ♦  m  X  -  ^*)|  •  -  J,  . 


(3«) 


•  ft 

atnc*  a  la  arolencth,  X  -  «*,  «dtcrw  K  la  tte  welt  nonal  vector  and  «  to  tte 
ewwUn  at  tte  ray,  felt Iplloat  ion  of  (52)  by  bjrltiU 


..m. 


A20 


r?«  •»  N*^(|  n*")  «  nW-'di, 


(33) 


We  nov  npply  (33)  to  the  incident  ray  at  the  turning  point  At  this  point, 
N  »  (1,0),  hence  ■  n*(xa)  und  K  -  n*  (x^/nfxg).  Thus  we  see  that  in  (51) 

K  is  the  curvature  of  the  ray  at  the  turning  point. 


fleyond  the  turning  point,  the  refracted  ray  is  given  by  (l6)  with 
xb  replaced  by  Similarly  the  phase  is  given  by  (17)  vith  replaced  by 
x^,  The  amplitude  on  the  refracted  ray  can  be  calculated  by  our  earlier 
method  although  a  technical  difficulty  arises  in  computing  from  the 
ray  fornula,  (An  Integration  by  parts  must  first  be  performed  since  straight* 
forward  differentiation  leads  to  an  indeterminate  form,).  The  details  vill  be 
omitted  here. 

We  new  consider  the  diffreotod  field  in  the  shadow  region*  lbs 
limiting  ray  is  tsneent  to  the  boundary  at  the  point  •  (xfe,  y'b).  It 
gives  rise  to  a  surface  ray  vhieh  procoodo  along  tha  boundary  x  ■  x^  as  a 
straight  line*  At  each  point  ■»  (x^,  y^)  the  eurfw^e  ray  shads  s  diffract  ad 
ray.  For  thsec  rayc,  a  -  n(xQ)  sin  ^  ■  r.(x^).  Honcc  from  (b)  the  diffracted 
rays  art  given  by 


y  -  y4(*) 


\  c»- « 


'ZFr 


(>*i 


Thus  they  form  r  one-parameter  family  of  oongruset  curves. 

In  order  to  apply  the  formula  (18.1b)  for  the  diffrueted  field  we 

must  evaluats  tha  limit  (l8*i;>).  From  (lb)  vs  sat  that  dy  •  ey. ,  hsne# 

0 


iv  ■  dy  cos  y  ■  dy^  cos  y,  Let  P*  ■  (x* ,y' )  be  a  point  on  the  diffracted  ray 


joiiJ.ne  P^  snd  P  ■  (x,y»,  Then 


da(P')  dv(P')  cos>(P‘) 

"  TVm  “  S'os  yfrj  * 


In  (18,15)  oQ  denotes  the  distance  froo  P^  to  P'  along  the  ray.  Since  7  *  */<; 


when  0  ■  0  ve  see  that 
o 


11m  cc«  y  l?' )  b  Urn  a0*  r^°)  _  d  cos  rl 

»t_*  “  J  <T  |o_  - 


P'-i^  o 


dr  „ 

•*3o0  • 


Mere  *  Is  the  curvature  of  the  diffracted  ray  (hence  the  curvature  of  the  Halting 
raw)  at  the  boundary.  How  fros  (18),  n(P)  eo*  r(P)  •  (n8-*8)1^- 


[n8(x)  -  n8^)]^8,  tbarefore  froti  (1B.15),  (35),  and  (36) 


tf(P,> 


cos  r(p|  1 


“M‘lf  -  lia - wa-m-fc-  -  — s— 

n(l»}k{Pj  n(P)»_  ®oa  r(P)  n(P)  eoa  r(?) 


*[»*(*)  -  n8(«^)j 


(10,11)  uow  yields 


ud(p)  •  sJdSj)  09  ^^(Qj)  ♦  f  ndo  »  I  b(v  ^ 
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dw^) 


Since  the  surface  rays  are  straight  liras  j-.j.t--  .  1.  Furthermore 

h  {l* 


n^)  ■  nO^)  and  ^  ndo  ■  n(s^)(yb-  y^' ),  The  surface  S  Is  a  plane  on 

which  n  and  7  are  constant,  therefore  the  diffraction  coefficients  at  P^ 
and  ^  are  equal  and  the  decay  exponent  ot^  is  constant.  Thus 

A 

*  OjtoJdo  .  (^.4)^, 

1 

To  calculate  the  change  in  phase  along  the  diffracted  ray  ve  use 
(3k)  to  Obtain 


d»  «  [dy2  ♦  dx^8 


— -  n  .  dx 


(») 


and 


(M'(  'l [A”^l]l/s “ *■( r^F5 


CnVf^n^8  dx  +  nU^Ky-y^. 


(<*» 


It  obtain  \(\)  «•  mat  uee  (31)  rather  than  (la)  which  la  lntefeandnate  at 

V  *■»  (kl) 

*J(\)  -  *lt/k  ^a(a«k)*^8  [«*(^>  - 


..... 

i* 

1 

r\ 
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Free  (9)  xo 

/(C^)  -  n(x^)y^  +  ^jf  [n^n2^)!^6. 

(U&) 

t  • 

We  say  now  insert  all  theoe  result#  in  (38),  This  yields 

to) 

ud(P)  **  QaS(x0)— **^(3^ [n£(x).n2(:tb)]’A/4  h 

/ 

a  *  , 

*  * 

* 

XX 

x  exp  /ihjn^JSr  **] 

t  * 

. 

The  diffraction  eoeffleient  e,  and  the  decay  expcaent  a}  are  obtained  frea 

(19,13  -  ID  by  asm*  •"Hx)  -  %  *<x)  -  *  •  «•<%)/«(*»,)»  “*  »(*>  - 

• 

Ihtti 

3  .  e*1^6  <ij  [(rHnf^)]^3  *^3, 

(«*) 

V 

end 

* 

to) 

- 

to) 

<*' 

the  profit*  discussed  lr  this  section  can  be  solved  exactly  by  separation 
of  variables.  The  asymptotic  expansion  of  the  solution  can  be  obtained  by 
•pplyio*  osyaptotlJ  aethni*  (the  nethod-)  for  orderln*  differential 

^lul.  \L  m,  »  «.  ««  ««  '*"•  -  *«  *"«- 

* 

r 

oo_ 


31 


rK 


B.  Aavwptctlc  ifrthods  for  Harwell 1 3  icnMtions 
B'l,  TlMB.tor»Mle  Solutions  of  tfcxwell'a  equations 

Although  solutions  of  the  rave  equation  ere  frequently  used  to  describe 
optical  phenomena,  It  is  veil  known  ihafc  a  rigorous  description  of  optical 
and  other  elect roengnetic  phenomena  can  be  obtained  only  by  solving  Maxwell's 

•mm  fry*  nrf  vr—nwaH r  If  MffK  fw^atiawrUt .  »«ywrif  e&if 

astluds  ore  particularly  useful  far  this  purpose.  We  shall  see  that  nany 
features  of  tbs  a.ytqptotlc  set  hod  for  salving  Maxwell's  equations  are  slnllAr 
to  those  which  we  have  exsained  for  the  reduced  rave  equation,  and  we  shall 
asks  full  use  of  the  similarity,  nevertheless  tbs  vector  character  of  the 
- 1  w  trees  gnat  ic  problaa  Introduces  significant  differences  which  we  shall 
— in  detail,  for  the  ttterlal  in  this  chapter  we  ere  largely  Indebted  to 
R.X.  Uawiburc  [3*)« 

Proceeding  ••  In  chapter  A  we  shall  a  muss  h*  nmolo  tlaa  dependence 
and  derive  tbs  tlae  reduced  fora  of  Itamll's  equations.  Then  an  asymptotic 
■cries  will  be  inserted  Into  the*  equations,  and  equations  far  phase  and 
mmiltuds  functions  will  be  derived.  Ms  will  see  that  the  phase  function 
sfpin  satisfies  the  decani  aquation,  fur  this  reason  such  of  car  earlier  work 
will  be  directly  applicable.  Zn  particular  we  stall  have  tha  aaaa  ray  aysstions, 
The  aula  difference  la  the  electro* gastic  tiwery  Ust  la  tha  vector  okwetw 
of  the  a^Utud*.  Bwewer,  even  here  the  essential  feature  renalnei  the 
eosponsnis  of  '■'<*  aapUtuds  functions  satisfy  ardiaary  differential  equations 
(transport  sqastlons}  along  tbs  rays. 


r 
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In  M.K.S,  units  [36],  fexueU's  equations  teXe  the  font 


v  (•£)  * 

(1) 

/S  \  A 

n  *£  ♦  £<»#)  -0, 

(^> 

A 

r»  4  4/  t  m 

v  •  u*  /r  ;  -# 

(3) 

A 

V.  (tg ) .  0  . 

(U) 

^  A 

Here  £(X,t),  7f(X,t),  are  the  (reel)  electric  end  Baguette  field  rectors, 
end  «(X),  u(X),  end  o^X)  ere  the  dielectric  "constant",  engnetic  pereeeMllty, 
•ad  conductivity  of  the  •edi'S'.  e(X,t)  le  the  electric  otsr^  density. 

*,  it,  cad  are  eisunsd  to  be  piece* rise  seooth  'unctions  of  X. 

We  shell  be  Interested  in  tlee.hereonle  fields,  of  the  fora 

f(X,t)  .  Rsj^  (Xje’^,  •  he[#(X)e-1,,t]  •  (5) 

Then  It  le  easy  to  see  that  (1)  end  (2)  ere  sstlcfisd,  provided  the  (eeeplev) 
rectors  satisfy  the  tine  reduced  equetlcns 

VxV-eins  £  •  •*£  I  (6) 

free  the  eeoood. equation  it  follows  lmdletely  that  v.(„  >f)  •  0,  so  {))  is 
euteustically  satisfied,  («)  any  be  thought  of  siuply  ee  e  definition  of  A, 

■Xb  our  aetetloo  w  hove  reeerved  the  ejnhole  t  end  8  for  the  lending  tens  of  the 
enplltudes  of  the  electric  and  aigrette  field  vectors,  this  eeeouate  ter  the 
uaorthodoe  notation  In  the  first  fee  equations. 


Here  A,  B  are  vector*  with  real  c opponents  end  the  her  ‘’-ruotes  the  eorrlex 
conju^pte,  Thuc  free  (5) 

/V 

C  ..  I  /  X  -  A  aa«  #.*  4.  »  «<n  l.* 

*  C  '  ”  —  . * 

Xt  follows  froa  this  equation  that  a*  t  varies  (at  oech  point  X)  the 
vector  ft  (X,t)  ■  (  ft  ft  2»  &  Ascribe*  an  #111  pee  vhich  Ilea  In  tiu  plane 
dr:eralnod  by  »iv'  *.  Tills  plane  of  polarisation  xa  therefore  perpendicular 
to  tia  vector 

BXA.^-  QBxA.  iAxg  .  JjgKt.  (9) 

The  principal  axe*  of  the  el  11  pee  correspond  to  the  e/treoe  values 
of 

ft*  .  $  (  V  t*"1*  ♦!*  e*1*  ♦  »  ft  •  fi  ).  (10) 

■quatlig  to  aero  the  t  Aerlvatlve  of  (10)  yield* 

.  ♦  Jv/L**  (u) 

-s 

If  we  Insert  (11)  In  (la)  vc  aee  thfct  the  extroao  vcluee  of  ft  are 

yn5"  >. 


(») 


f 


0 ' 


VV’ 


"  5  («  E  *  *H~) 


and  the  encrpj-  flux  vector  (Poyrrtisi;  vector)  Is  defined 


A  A  A 

s  .  Ex  V. 


He  oenne  toe  corresponding  tiw.averae**  quenlilivn 


v £  («£  2  ♦ 


■■«  / 


£  ^ 

(£  xtf  )dt. 


Drat  (8)  ud  tbe  awlogra  equation  for  "U  it  ia  way  to  show  that 


Iml* 


<■  4 


j  u  i  yoaltiva  lot  ajar,  or 


t  ee, 


Zb  wvty  ayoM  «(X)  end  p(X)  have  the  oenateak  ealnaa  t„  •  8.8jkX10r' 
AuMd/natar  end  uQ  -  I.STIXIO’*  henry/ —ter  ()ti).  tbe  a  patent  o0  •  (t^»0)*^8 
•  S.99T9O0f/°  Mar/aaa  la  the  taalUar  *ayoad  of  Ugbk*.  Aa  ia  obaptar  A  va 
lakro*MM  Uta  yropapMaa  eonatonk  (er  *i«a  aatar)  k  »  a /e^,  and  w  a  awe  ttrt 


» »*i #.* *vth- S'  f?  S’ 
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re 


the  complex  vectors  <?  ,  ?f  have  ssysptotjc  expansions  of  the  fom 

y  (no-  ^  nr*  v  (i) 

*55  Bs»o 

The  reel  scalar  function  s(X)  is  again  called  the  phase  function  (or  phase). 
If  vt  Insert  (l)  In  (1.6)  and  collect  coefficients  of  the  sane  powers  of  (lit) 
we  obtain 

* » K* "  *  Vi  •  w 

(?) 

•  •  Mum  • 

The  equations  for  ■  ■  0  are 

8  ♦  o  il  ■  3,  S’*  Ml  •  o  til  ■  0. 
o  o 

•are,  and  la  all  aufc  sequent  aquations,  we  unit  the  subscript  sera. 

Ms  sea  ait  ones  free  (3)  that 

M*l*«*ft^*0  . 

hr  ellndnetlng  I  frca  (3)  we  obtain 

C02t|l  *e-%X(^»yf)w  {%)*«.  (j) 

It  follow*  that ,  if  S  ;•  ie«t**«ro,  *(x)  oust  atf:^ 

the  alooaal  squatloe 


.  99  - 


v  •  jf  I  E‘l  ■  W  H’H, 


S  -  <*!)*  ,  -  -S-  9s, 

o  on 


ee*  »  5  •  **, 


In  order  to  obtain  differential  equations  for  I  end  B  elans  •  rey, 
™r  rotam  to  (C.2)  for  ■  •  1.  Tw  convenience,  we  set  'fc  ■  1  and  ayraatrlae 
the  equations  by  introducing  e  fictitious  negnetlc  conductivity  tj  (which 
latter  will  be  act  equal  to  aero).  Than  (2.3)  yield* 


*1  *  V  »  e^‘V  -v  **  -  •afB* 


etloea  era  eyanetrlo 


^  -Uj,  -*j. 


Ttoey  cine  Ugly  new  uoulitisaa  on  *  and  It,  To  Obtain  theaa  conditions,  m 
first  note  that  K  k  (R  x  g^)  •  (►tl)»»j®e**1  #  alma  ^  •  «J*».  Xf  «e 
mltlyly  the  first  aquation  (9)  toy  e^  aad  add  1%  to  the  vector  product  of 
with  the  i wood  squat lae  «e  atotale 


(fcS^)l  *  -  Xl.lxJVxll.yxS'Y  Sj*.  (10) 


'  IOC 


B3 


But  trm  (2.3) 


R  x  B  ♦  cetE  ■  0,  SX1.  c0uH  -  0, 


(U) 


Brace 


...  r 


n..  #1  •  ,  1  «  m  ..  «» 


,V 


J 


(12) 


VT  T  V“j  ’  v* 


Tbua  w  Lara  oUUlmd  on  option  Involving  X  froa  vtaleh  K^,  B^(  end  B  1.  ;ra 
beuii  •Xlndmted, 


Bow  V  x  X  -  7x  (V  ■)  .  0  so  Vx  (i  R)  .  7(i)  x  I, 


IX 


(13) 

^  tx  (Vxl  R)J.XX  | *x  j^(i)x  *J  j  .IX  |  -|r*ti)]R|  - 


(IS)  ray  *•  written  u 


(U) 


■  *x  (3  I  x  X)  ♦»  X  x  f*x  I)  ♦  I  *  <v*  j;  *)  -  «e(^  ♦  #^)*j  •  0, 


Mt  BOV  Mt 


1 1  •  A  tad  «m  tte  raster  lkttty 


(13) 


9x  (A  X  I)  ♦»x(Vxl)tl*(9x»).  -H\  ^  ♦  Rj  ^  ♦  Aj  ^) 

.  X*A  ♦  tf&S  t  V(A«8). 

►t  *  o  it  follow  trat  (it)  ts 


•  Of-jUkit**. 
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(16) 

+uT>(iR)  4e0(i#J+^)lj  -0. 

let  X  ■  X(t)  b«  the  aquation  of  u  ray,  and  lat  ua  chooaa  tha  paraaatar 
t  ao  that  |X|  •  |j-|  -  a.  [Sea  tha  ray  aquatlana  (A2*8),  (A2.9)j.  than 

Xm%mK.  (1?) 

thus  (R*9)B  ■  (X*V)£  »  ~  ,  Ua  um  thla  In  (16)  and  not*  that  tha  quantity 
In  bnar.aa  wet  ba  yaralla'  to  R.  ffcarafora  (16)  aay  na  vrlttan  aa 

H  ♦  Jj  B  R)  ♦  -^(t  «8  ♦  Ma^)*  •  fli  •  (Ifl) 

Sara  P  *»  to  -a  Ictcrmaod.  awarer,***  *  n*  a***  mg  •  0,  Tharaftwa,by 
aoalar  wit  lyll  cation  of  (IB)  by  *,  wa  obtain 

*-i«  (*•«)-?  (*•!)•  <!9> 

Mrthotftra,  turn  (AM) 

(•*)  .a  (80) 


(«) 


Aatloa 


Ho*  (19)  aad  (10)  ytaU 


-  ,'.02  - 


u*(i  R)  -  (J  *  Xj)Xj  -  V  <22) 

(if  m  «•  const,  Aj  Is  the  Le.placi*n  operator),  How  (18)  becomes 

J?  +  |  33  V  +  (*?)  7e  +  5s  (,52  +  W01)E  "  °*  (83) 

By  swans  of  the  symmetry  property  (9)  we  obtain  on  equation  analogous  to 
(23)  for  H.  Vie  then  set  o£  -  0  in  (23)  and  in  this  equation  to  obtain 

s  *  1 1 V 4  sv1* 4  (¥*)  <“> 

S4i<V4w>"4  ft?)  «> 

'xbese  ordinary  differential  equations  for  *  and  H  alon*  a  ray  can 
be  simplified.  7c  do  bo  vc  introduce  the  vector 

e  2 

,-«uR-^  "*•  (aS) 

Than,  elnoe  $ylag«»|?«»X.»|?c»Ra»tr»V«# 

2  2 

ll  >  ifr(tf)  •  3*  ♦  MP»^»  ■  fig  ♦  $—  log  «  .  (a?) 

u  o:  c” 


y 

'ft' 


/ 
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fly  inserting  (27)  into  (24)  ve  obtain 


thus  (26)  becoma 

±  (  V«  *)  ♦  $  (jjj  ♦  o^n)  /«  B  ♦  *  •  0.  (30) 

The  amln»nw  equetiea  for  II  ia 

J;(/5»)  ♦|(jSr*e^u)  Vb.o,  (31) 

9mm  equations,  otOcb  dftmdnt  h <*»  t  tad  I  vary  alone  a  ny,  oon 
bo  wylaaad  by  oUflUr  equations  far  tbo  anyiltudo  and  direction  of  th»M 
tartar*,  9s  this  end  m  introduce  a  real  Malar  function  defined  by  tbo 
differential  aquation 

lc«  »x  -  (4j  V  •  f  ♦  e^u)  (J2) 


and  U»  Initial  condition 


-  loU- 


B3 


2v^  m  2v  »  «E*1  ■  it  t  »  to  .  (33) 

Here  tq  la  a  one  point  on  the  ray.  Lei  r,  Q  coqple*  veclora  <l«riiMu  uy  the 

equation* 


/«  S  »  /S^  P,  yS  H  •  /5Ta  Q  %  (3*0 

Than  P  and  Q  arc  unitary  vactora  (!,»,,  P.f  •  Q*S  ■  1)  to*  t  ■  r0.  Whan 
(3*t)  ia  Insert  ad  into  it,  (30)  baeoaaa 

(35) 


1 

3 


/d  log 


and  (31)  la* da  to  a  siailar  aquation  for  Q,  Thia  aquation  and  (39)  oi^Llfy 
•hoc  (38)  ia  ueed,  and  th*  roaulta  arc 

(36) 


m  0, 


Wtm  (3*0  and  (2»  ve  sea  that 

P*Q  •  P«*  •  **•<!.  (37) 

If  ee  Multiply  the  firet  aquation  la  (36)  by  ft  and  nan  the  feet  the*  ?.  •  0, 

ve  detain  ?  •  ^  *  0,  Us,,  fc?  •  eoMt,  In  the  sane  vay  ue  prtwe  that 
5*Q  •  oooat,,  and  it  follows  that  P  and  Q  are  unitary  veotoye  far  all  t.  It 
alee  foUoee  me  (3*)  that  an»i  •  p»F  •  Than  me  (>)  «e  eoe  that 


B3 


* 
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v,  •  v  and  «t  my  hanc (forth  ocf,t  tba  «ub script  1«  fha  differential  equstions 
(3 6)  vill  ba  further  analysed  in  section  6,  Mar  ve  shall  easaine  (32)  vhlsh 
datendoes  the  saro- order  average  energy  density  v. 


Since  log  v  ^.Vs,  (32)  yields 
♦  vJai  ♦  ♦  ee  o^ilj  •  0 


^(■^  ^»)  ♦"5  o^u  ■  0* 


Me  now  set 


«■«« 


and  note  that 


.  £  -  e^u  q. 


(»)  *nd  (hi)  yield 


*(d  ^  *>  "  «(*  <J|  Vi*]  "  *  <*> 

But  (h2)  la  of  the  wane  fans  cc  equation  (Aj.h)  and  haaaa  n  slatle  agpdlaatlcn 
nf  (koae*  iheor— ,  as  In  MQtler  \J,  yield* 


d  . 

a 


I* 


r 


r  ai'Sitir •*-■*'  *  -Ji 


-  icA  - 


If  o  la  an  arc  length  para  Deter  on  the  ray,  (At>*lu)  shows  that  do  ■  ndr.  Hence 


e  Md®  rr 

e0ndr  ■  -*■ —  do,  and  (>A)  becomes 


“38“ 


»(»0H(ae) 


{•  (  ^  -} 


Equation  (43)  Aeternlnei  the  variation  of  the  tero- order  (vtiagt 
energy  density  v  along  a  ray.  It  Is  the  analogue  of  the  solution  (A3. 8) 
of  the  xero- order  transport  equation  for  the  rednoed  saw  equation, 

((e)  »  |||^|  Is  the  expansion  ratio  introduced  io  section  A3,  The  hlgliev  wier 
transport  equations  for  IuonU's  equations  are  anUy  «d 

Since  »  •  »^  to*  values  of  the  aero  order  field  vectors,  I  an*  q  ean 
he  deteralnad  from  (3*0  end  (1*5)  ones  the  polarisation  vectors.  P  and  Q  are 
found,  the  equations  far  P  and  Q  ere  studied  la  the  next  section, 

Xa  a  nedluw  for  which  e^  a  0,  (k?)  beconas 


ooneervatlon  in  a 


due  to  the  oonduotlvlty  of 


Mi  equation  enreeees  the  utUUkaorn  oriuolple  of 
td he  of  raye.  (kj)  describee  the  dissipation  of  aaari 


the  aadlua. 


According  to  (1,9)  the  plane  of  polertaatioa  sa 
vaster  $j>£x  i ,  so  to  taro  erdar  ftt  la  perpendicular  to  1 


to  the 


cular  to  the  vector* 


9 


-  1C? 


BU 


|jExIand|jPx’P,  But 

X  (P  X  P)  .  (?a  •  P )i>  -  (^s  •  ?)P  -  0.  (1) 

Therefore , the  plane  of  polarlaatlon  1»  perpendicular  to  perpendicular 

to  the  ray,  iron  (1.13)  ^  m«  that  to  taro  order  the  elllpticlty  la  given  by 

j,  u-VR*  i- 

'  ‘  »i .JFT  '  i.JF?  '  (a) 

■guetlon  (3.36)  lapllea  that  P2  and  ^  are  eonetant  on  a  ray.  Banoe  the 

aUtetloity  ia  content  airy  «  ray. 

Wtm  (1.17)  «•  Me  that  for  t&UULfit i tMSfc JBiaflgttai  f  u 
proportional  to  a  real  vaster*  i.e.  * 

(3) 

then  ?0  la  real*  and  a  la  a  ravlax  water  of  ssdulue  we  vhleh  la  eonetaok 
on  a  ray.  l\otberacre  I 

1  -  9  .  f  -  l§  ,  (t) 

l.e. ,  P0  la  a  ml  adt  vector. 

Proa  (t.3)  and  (3.34)  ve  new  eae  that 

Q  .  t  *  P  -  •  V  (5) 


* 


vhi»rc 


“  10ft 


7  -  -s-  ^  \  <*..  (6) 

1*1  “ 

T,  r  .  and  Q  oro  orthogonal  unit  root  or  s. 

o  o 

ftrtherwore.  it  is  ensy  to  sea  that  to  zero  order 


M 

cosfVa-wtjP^ 

(7) 

£  fz 

ra  -e  , 

jri~i 

cwJjKS-VTjj 

V-/ 

Sere  ve  have  absorbed  the  (constant)  phase  of  a  into  a  which  la  undat trained 
up  to  an  additive  constant  on  a  ray, 

Tor  the  oaaa  of  linear  polarization  ve  any  replace  P  by  PQ  lu  (3.3^)> 
and  write  that  equation  in  the  fora 

BP’  ♦  P»0  •  («X,),JX’  *  °*  (9> 

Sara  ve  have  need  (A2. ')  In  the  font 

(•*’)•  .7a.  (10) 

the  prim  denotes  dlffartHtlatlon  with  respect  to  the  arc  length, fl.  However, 

Vq'  X1  •  0.  Hence  F0*  (nX* )'  ■  P^fnX*  *  n*x')  *  n^o*x"»  4,1,4  beemea 

?o  4  <V*")r  *  °*  tU) 
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Bit 

We  next  apply  the  theory  of  space  curve*  to  the  ray,  und  Introduce 

X" 

the  taneent  vector  T  »  X',  the  principal  noral  vector  N  ■  —  ,  and  the 

lx") 

binoraal  vector  B  -  T  x  .N.  These  vectors  satisfy  the  Ftenet  equations 

T»  -  <N,  (12) 

H*  •  .Kg  +  yB,  (13) 

a*  ■  -r«. 

Bare  *  ia  the  principal  curvature  and  y  la  the  torsi  or.  of  the  curve.  With 
these  foraUna,  (U)  becomes 

P&  ♦  <(P0  S)T  .  0.  (13) 

Since  PQ  Is  nonal  to  T» 

P0  -  *  ♦  SB|  0?  ♦  0* . 1.  (16) 

If  ue  Insert  (15)  Into  (13)  ve  obtain 

0PB  ♦  OB'  ♦  B'B  ♦  SB'  ♦  «B  -  0,  (17) 

or 

(<*  .  rt)B  ♦  (S'  ♦  oy)B  •  0.  (IB) 

Xt  follow*  that 


t 


a»  .  ys  »  Oj  S'  e  ay  •  Of 


(W) 


110  - 


4* 

s 

3 


or 


4;  (a  +  IP)  +  ir(c»  +  ip)  -  *\ 

UV 


a 

•if  ydo* 


This  equation  has  the  solution 

a  IP  -  (o^  ♦  130)  e  "°o 
Let  0  be  the  angle  between  ?Q  6*1  H  (see  figure). 


Then  (21)  becoaes 


and  from  (1 C) 


0.04 

O 


P  .  a  coe  0  •  B  sin  0. 
O 


(?o) 


(21) 


(22) 


(23) 


Since  Qq  ■  1  *  pQ,  T  x  B  «  »,  and  7  x  *  •  «»,  it  follows  tram  (2$  that 

Qg  *  I  sin  Q  41  eoe  0.  (24) 

(22).  (23).  anA  (2>i)  give  the  rrtt ton  of  the  ?ol*ri«et1"n  wtoes  T\, 


4 

5 


I 
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ralfetiva  to  B  and  W.  If  the  ray  reoelns  In  one  plane,  then  7  ■  0,  and  6  la 
coneWal  «iuu*  a  ray.  A  sufficient  (hut  nr*  *  necMury)  condition  for  this  la 
that  the  Medium  be  hanogeneous,  1. a., that  n  ■  constant. 


B5  Reflection  end  tranaalialop  at  an  interface. 

In  this  section,  ve  focus  our  attention  on  an  interface  or  surface  8 
vtdch  separates  two  redone  In  which  c  and  u  are  snocth  functions,  these 
functions  my  hava  Jus®  discontinuities  across  8.  In  raglom  1  and  8,  we 


The  values  of  the  refleeted  and  tranatftted  fields  et  the  interfere 

oen  be  derived  free  the  vell-toown  di scontinuity  conditions  far  the 

elect rorngimic  lisiu,  •UiuU  twtul.o  the  essiissity  tf  mwwM 

A  A 

of  £  end  #  end  therefore  of  &  and  Jp  .  Since  we  assure  that  the  lneident 

end  inflected  fields  ere  defined  in  rectos  1,  end  the  tmaMtled  field  ie  defined 
In  region  2,  the  conditions  beecne  (l) 

»x<£i ♦  *x<v1^>*»x^’t*®#* 


4 

if 


% 

* 


i 

« 

* 

4 

a 

* 


tore  a  denotes  e  volt  wester  none!  to  r  poAatiac  in  the  direction  tree 
red  on  1 1®  redo®  8  (see  fi«ure). 


-IIP-  H5 

B»ch  field  is  of  the  form  (2.1),  so  it  satisfies  (2.2),  (2.3),  etc. 
By  inserting  the  expansions  (2.1)  into  (1)  ve  Obtain 

a.  ■  s  »  s. ,  on  fl,  (2) 

1  r  v 

nn»l.  Tfvr  t.h*  <*o*fftrrt«n[tBs 

H  x  t\  *  Bp)  •  »  *  V 

Rroi  (2.3),  we  haw 

»1»VVA"' 

WA-0' 

via  now  introduce  the  paxaactric  equation  Jar  the  wrfhoe  8, 

x .  x(^,  -  x«).  (?) 

then  (2)  my  be  written 

-  S0W3-  w 

DtfferentlrKao  of  (8)  with  respect  to  ^  sad  (g  yields 

’  *t  *  *r  *  ^ 

mm e  the  waters  X,  ere  ta^saMsl  te  8,  (9)  i mfiim  that  the  differeoeee 

*1 

^  ere  norsel  te  8,  l«e., 


8  J<  ♦  8*.)  •  »  *  V  00  8‘  '3J 

"i”*!- VA'0*  <*> 

%  *  «r  '  VA  *  °*  °> 

\  -  «0I*A  •  C*  (6) 
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\  *  7t*,  *»t  -  \  ♦  (10) 

It  follow »  that  ^»r  and  (and  hence  the  reflected  end  transmitted  rays) 
■»t  8  lie  In  the  plane  of  Incidence  determined  by  and  R.  (Tils  plane 
la  the  plane  of  the  figure. )  Furthermore  (10)  lilies  that  there  exist*  a 
unit  vector  V  and  a  real  scalar  a  such  that 

^xi.^xi.^xi.al  (11) 

V  la  perpendicular  to  the  plane  of  Incidence,  (in  th->  figure  V  palate  into 
the  paper). 

Sy  equating  the  aagnltudes  of  the  vectors  in  (11),  we  obtain 


n^  aln  ■  n^  aln  0  •  Dj  ala  O'  •  a.  (12) 

The  angle*  0,  O',  0*  appear  In  the  rigur*.  Xt  la  clear  that  the  only  eolations 
of  (12)  consistent  with  the  figure  ere 

0"  •  «  -  0,  (13) 


aln  0* 


da  0|  0  <  0*f  g/?. 

k 


(u) 


(u)  and  (It)  say  be  rcaci^lsd  ::  t*^  1>“  of  rtf!*-****,  and  RaH1!  ton 


J 

i 

i 

e 

•i 

\ 


i 

* 

3 
» 

S 

*■ 

■5 

4 

i 


^  da  0  >1 
"8 


(15) 


then  (lb)  has  no  real  solution  O'.  This  is  the  cose  of  total  reflection 
discussed  in  Section  A13. 

In  order  to  detsndns  the  amplitudes  of  the  reflected  and  t  re  emitted 
fields,  it  is  convenient  to  introduce  the  3  unit  vectors 

Vi  -  V  x  Vs^,  V.  -  V  x  Vt  -  V  x  (16) 

which  appear  in  the  figure.  81nce  is  artho*.oael  to  It  oan  he  expressed 
as  a  linear  ccnhlnstion  of  and  v.  The  setts  aeaertlon  applies  to 
1^,  t*,  Hy,  1 Thus  we  my  set 

/«A  *  Vi  *  aiT»  A*r  '  'V  Vr  *  ®rV»  -  Vt  *  V*  (l7> 

Than  It  Is  easy  to  shew  that  (b-6)  are  set  !*<*!•*  je'M'ied 

M 

^1,  ■  •»,*,  ♦  v*  ^  *,  •  •  W  V>  1%'AW' 

Xu  vte  to  apply  (3)  vs  first  Observe  that  V^x  I  ■ 

(f  X  ^/b^)  X  S  •  -?0%,  •  V1^ )  *  •?  0,  and  therefor* 

-  Q^M)  ♦  .  -O.f  cos  0  ♦  Pt««.  (19) 


-Ilf  - 


25 


Similarly. 


(a) 

<x  ».  P_  a 

x  H  ■  — S—  V  cot  0  ♦  —i —  V  x  B,  Hrx»  -  .  -s—  V  cot  0  ♦  — £ —  ▼  x  », 

^*1  ^*1 


-<*.  *t  \  \ 
Lx). V  cot  O'  ♦  — i  V  x  V,  *.x»  •  — i  Te«0'  ♦-»TX«. 

^  t  ^Ma 


Bart  ve  !*v«  uMd  tht  Met  that 


eac  0“  •  -r(x  -  0)  •  -cot  0. 
If  x*  bn  lM«t  (20*22)  loto  (3)  «•  obtain 


#r**l**^j  Sf 


(23) 


(2k) 


r  -S 


4 

$ 


^1 


3s# 


$ 


/jl  cos  0* 

J  Eg  COS  6 


f*l  J.  COB  0( 

V  ^  V  **2  COJ  °" 


f25-B) 


The  conponents  a,  g  of  the  fields  are,  respectively,  parallel  and 

..  1  1.  4U*  4«W4mma  ▼+  4  m  A,%MA44«M>e  Axe^nwnan  4a  »>r*  fk* 

UW4UW*  WV  yttu  V*  *«  »“  j  •  ■  ""  *  •■*■ 

notation  for  the  parallel  and  nomnl  conponents  of  the  electric  fieldi 


a.  a 

1  I* 

E.  *  .  il  *  * 

«  y-  •»  yq 


If  ve  assume  t.tict 


'Btp«  i'Bin‘ 


»  s™,  *  ■„“» 


■  u.  «  u 


Uj  ■  u- 


ond  set  n,  -  cc>/tp  ,  u,  *  c=  /c.u,  then  (??)  leoowes 


cos  0* 
^  '  n8  cos  Si 
Eiy  ’  Eip  "  cos  0* 
1  4  n'2  co  s  Z 


m  a  V 

Hp  "  "<p  n^eo-i  u*  * 

^  +  rTeoI  0 


°1  coe  01 
“  cos  t> 

E«»  *  ■h  o»  * 

4  eo«  A 


V  ■  Ilt,  t-8 - . 

tP  lB  "l  .  co.  O' 
£  4eoT5 


V 

U 

N. 


IV 


t 


-  117 


B6 


i 


These  am  identical  to  the  Fresnel  formulas  for  reflection  end 

transudes' on  of  a  plane  electromagnetic  wave  at  •  plane  interface.  We  have, 
of  course,  shown  that  they  ere  valid  for  the  sere- order  tame  of  the 
nsyqptotlc  expansion  of  an  erbltary  electromagnetic  wave  at  an  arbitrary 
(smooth)  interlace.  ay  using  i'uw  r«»ul«.»  (C“)  ns  initial  condition*  »•*«■ 
electric  field  on  the  reflected  and  transmitted  vays,  the  tero-order  reflected 
and  tranaaitted  field#  can  he  found  away  froa  the  Ibterlboe. 

B  6  Reflection  from  a  Perfectly  conducting  eurface 

The  well-known  condition  on  the  elaetrcngnatlo  field  at  the  surface 
of  a  per  fact  conductor  ie  that  the  tangrrtlai  component  of  &  wet  vanish. 

Zn  contrast  with  section  5,  we  have  only  Incident  end  reflected  fields,  end 
the  boundary  condition  any  he  stated  in  the  fern 

*  *  (*i •*  V  *  w 

The  consequences  of  this  condition  can  be  Obtained  easily  by  eUply  nodi  tying 
the  equations  of  section  5  in  an  Obvious  wey.  Zn  this  section  we  list  the 
nodi  fled  equations,  using  the  seas  equation  maters  to  ;*cilitate  the 
•aagmeon.  thus  we  Obtain 


*t '  V  0,1  a* 

(7) 

V  *  .  0 

(3) 

0*  -  «  .  0 

(13) 

i 
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i 

Ei "  aivi  ♦  piv> 

yft  Er  .  apv,  +  Prv 

(17) 

; 

v/S  \  -  +  ^V, 

*/5  Rr  «  -0rvr  +  ay 

(r8) 

j 

1 

\ 

°r  *  “i  * 

0r  *  -Pt, 

(35) 

\ 

4 

i 

"  8ip  * 

*rn  "  "®in* 

(28) 

Bare,  aa  in  section  5>  subscripts  p  and  n  denote  component*  of  end  Ep 
parallel  end  nornel  to  the  plane  of  incidence.  As  before*  these  seines  can 
he  used  as  initial  conditions  to  determine  the  reflected  field  all  alone 
the  reflected  rays. 

87.  Ihdlstion  frost  sources,  diffro**-*"*,  f»wy 

la  order  to  discuss  radlstion  fro*  sources  sad  diffraction*  only 
slight  Modifications  of  the  results  nf  chapter  A  are  required.  In  this 
section  we  present  those  nodi  flections*  tocsthsr  with  s  sundry  of  the  results 
of  the  present  chapter .  Equations  vhioh  are  taken  froo  earlier  eeetlons  base 
augers  to  tnair  left  vhioh  Indlant*  their  arista. 

Xa  M.K.S.  units*  ILu  real  Ussr«lstr«onio  alectric  and  angnetio  fields 
are  clvsn  by 

(1-5)  i(X,t)  .  Rs[e  <«**<],  £(x,t)  .  B.JjftDe*1*].  (1) 

fbe  eoopltx  vectors  $  (X)  end  Jaf  (X)  netlsfy 


j 

i 


i 


a 

i 

« 

I 

> 


i 

1 

( 

i 

■i 


c 


(1.6)  ^xy+lw«£  »  ,  Vx5  rn  0.  (2) 

Here  «(X),  u(X),  and  o^(X)  or*  eiven  functions  vhich  characterise  the  nediun. 

8 

for  larce  '*  ■  u/eQ  (cQ  ■  3  X  10  meter/ sec)  we  Introduce  the 
asyndetic  expansions 

(s.i)  y  («>"«„,  <i»r*  v  (31 

ttz  fee 


The  tero  order  amplitude  vectors,  BQ  ■  E  end  Hq  •  H  satisfy 
(2.3)  ^XH4«a.O,  ^Xl.c  UH.O,  (k 


(2.t)  I-H  -  *.^i  .  g.%  •  0, 


vhT.?  s(X)  Hutlsfle*  the  eiconal  ecus t Ion 
(2.6)  (*)2  .  nZ(X). 


(?.?)  B*  ’  eo8,tt  “  v0 "  Ax)  1  c8(x) "  smSm  * 


Bquetloa  (k)  shows  that  I  oan  be  Obtained  from  a  knovladgs  of  B 
(and  vice- verse).  (5)  shove  that  S  sad  B  are  aituclly  orthosonal  and  «aeh 
orthocenal  to  i.o.,to  tha  ray.  fros  (6)  It  follow  that  the  phase  e(X) 


I./ 
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and  rays  are  given  by  the  equations  of  sections  A 2  and  A6. 

w»  define  the  zero.order  energy  density  function 

(3.5)  v  .  |<  E.g  r:  |u  H-I,  (d) 

and  the  polarization  vectors,  P  and  Q,  by 

(3.3»0  E  ■  P,  H  -v/f*  Q.  (9) 


Then 


P-P  •  1,  Q>Q  -  1,  P*2  -  0, 


(10) 


and  P  and  9  each  Batik*..  the  first  order  systen  of  ordinary  differential 
equations 

(3.36)  geS^L-^.o.  (U> 


Bare  o  denotes  arc  length  along  a  ray.  The  value  of  v  Along  a  ray  Is  given 


by 

(s.  *:) 


w(s,;  i(t ) 


Bare  (  (o)  -  Is  the  expansion  ratio  introduced  in  Section  A3,  Xf 

B  and  n  are  given  at  sow  point  oft  on  a  ray,  then  at  this  point  v,  P,  and.  Q 
osn  be  Obtained  free  (C)  end  (9).  A*  any  other  point  o,  on  the  ray  v  is 
given  by  (19)  and  P  and  Q  osn  be  Obtained  by  solving  (ll).  than  finally 
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BY 


at  a,  E  and  H  are  gives  by  (9). 

The  plena  of  polarisation  of  the  electric  field  ie  perpendicular  to 
the  my  and  tho  elliptieity  is  constant  on  a  ray.  For  the  special  ease  of 
linear  polarization,  i.e.,aero  elliptic  it  y,  additional  conelueioneoan  be 
drawn.  In  this  case, 


SI.  a  % 
PO/i 


t\.  f\ 


n  ..  aA 

*  "  “*0'  **  "  “■'o* 


\  * 


Here  a  is  a  coop  .ax  n  jitber  of  nodulun  one  aal  P?  and  are  real  unit  vectors 
voich  are  mutually  orthogonal  and  orthogoisl  to  the  ray.  Furthermore  to  zero 


order 

(4.7) 

£  ~  coa|)i3-vt]Prt, 

(14) 

(4.8) 

eo#pce-wt]Q0. 

(13) 

In  addition  (for  the  ease  cf  linear  polarization)  the  rotation  of  the  polaritttioo 
vectors  PQ  end  are  given  by 

(4.23)  1’0  *  K  cot  0  -  I  eln  0,  (1$) 

(4,  ?4)  ^  ■  H  ein  Q  ♦  B  coe  0.  (17) 

In,  |  and  B  are  the  unit  notrail  end  binorml  vectors  of  the  my,  and 

c 

(4.29)  0  .  Qe  ♦  /  Me’.  (18) 

*0 


s 


s 


r  ie  the  toreloc  of  the  ray. 


•*’22  -  - 

The  conditions  for  reflection  end  transmission  et  an  interface  are 

given  by  (2),  (13),  (it),  and  (25)  or  (26)  of  section  5.  Similarly  the 

condltiona  for  reflection  at  a  perfect  conductor  ere  given  by  (2),  (13'  and 

(25)  or  (26)  of  section  6. 

Homogeneous  media 

Let  ua  call  a  medium  homogeneous  if  n(X)  «•  constant.  Since  in 
applications  j*  in  almost  always  q  constant,  constancy  or  n  means  nnat  • 

is  constant  too.  In  this  case  the  rays  axe  •ti&it&t  lines  and  the  phase  is 

given  by 


(Ab.2) 


n  m  *_  ♦  r.o 
o 


(19) 


troa  (11)  va  see  that  -cue  vectors  P  and  Q  are  constant  on  a  "ay  and  hence 
(to  tero  order)  the  direction  of  the  on  jar  end  minor  oxen  of  the  ellipse 
of  polarisation  are  constant.  The  expansion  ratio  is  given  by 


(Ah.U) 

and  hence  (12)  becomes 


.  (V  °KP2  ♦  •) 
12 


(20) 

(a) 


(P,t  O  (p.t  c  )  f  r  r  ~) 

v(9)  “  v(oo)  ‘Ml  Vj  • 

and  P^  ore  constants  on  a  ray.  PTcw  (9)  and  (21)  v*  new  have 


B(«)  •  E(oo) 


!¥*$# 


1/2 


•xp 


^  -\A[  °i  d°'| 


and  a  similar  aquation  for  H, 
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Radiation  from  sources 

As  In  section  A?,  point,  line,  and  surface  sources  my  be 
characterized  by  giving  the  values  of  a,  E,  and  K  (or  of  s,  w,  P,  and  Q)  on 
the  source  manifold  M,  This  Is  particularly  conve:iient  when  M  is  a  secondary 
source  such  as  occurs  in  reflection,  transmission  and  diffraction.  The 


If  M  is  a  (non-characteristic)  surface  5,  then  on  every  outgoing  ray 
v»  my  measure  e  from  8  an.,  (c)  is  given  by  (23)  with  oq  replaced  by  0. 

If  M  is  a  point,  let  dO  be  an  element  of  solid  angle  or  the 
starting  directions  of  tbe  rays.  Then  far  »e-»  0  da(«0)  ~  ®0  d  Q  and 


7(0).  e?w(«).  (23) 

For  a  homogeneous  medium,  dn(c)  -  o*d fl.  Hence 

*(«)  ■  *0  ot  l(o)  -  (26) 

Here 

Can 


i 

l 

s 

i 

% 

a 

i 

1 

I 

I 

i 

4 

m 

• 

i 

* 

l 

i 

i 

\ 

I 

J 

i 

4 

* 

n 

3 

i 

1 

! 

1 
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« 
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Here  B(0)  is  ilven  by  (31). 

Diffraction  by  edges  and  vertices 

As  in  .action  Aik  if  on  electromagnetic  vave  (3)  is  Incident  upon 
an  ed^e  or  vertex  M,  that  nr. ui fold  acts  as  a  secondary  s<>urce  producing  a 
diffracted  vnve.  The  ^«ses  of  incident  and  diffracted  'raves  sulisiy 

a*  ■  e1  on  M,  (35) 

i4  -  (d)E3'  uu  M.  (36) 

The  aiffraction  coefficient  (d)  is  a  matrix.  As  in  section  Aik,  (35)  inpij.es 
the  law  of  edge  diffraction.  Tor  a  homogeneous  medium  the  field  diffracted 
toy  a  vertex  or  edge  is  given  toy  (26)  or  (30),  l(0)  toeing  given  toy  (36). 

for  an  lnhaaogeneous  medium  ve  ahnll  discuss  diffraction  toy  an  edge. 

The  discussion  for  n  vertex  is  similar.  Wo  first  usa  (29),  (6),  and  (31)  to  Obtain 

-d(o)  -  <"d(">  -  J;  -  |c(0)  E^O)-  ?(0).  (37) 

Then  Iron  (9)  va  Obtain 


How  on  the  diffracted  ray  v*(a)  la  given  by  (26),  (37),  and  (36),  while  F*(«) 
Is  determined  by  the  cysts*,  of  differential  equations  (11)  and  the  initial 


12b 


37 


condition  (78}.  Ha  vine  determined  vd(n)  «nc  ra(<»),E?*(c)  Is  given  by 


The  phase  of  the  diffracted  field  Is 

a 

od(c)  m  s i(0)  +  f  ndo,  (1*0) 

0 

and  the  field  associated  with  the  edge  diffracted  ray  is  given  by 

j,tryw^M.  w 

For  an  edge  of  u  perfectly  conducting  thin  screen,  the  diffraction 
coefficient  natrix  le  given  by  equation  (A12)  of  [£0]. 

* 

Diffraction  by  a  t Booth  object  in  a  honorentous  asdlu^. 

The  description  of  the  phase  functions  s°  and  ad  and  the  reyt  for 
both  the  surface  wave  and  the  diffracted  vave  is  identical  to  that  given  in 
ohapter  A.  In  order  to  describe  the  amplitude  vectors  on  a  surface  ray,  we 
first  introduce  the  vectors  2^,  the  unit  tangent  to  the  r»y,  Dj,  the  outward 
uult  »rwl  to  the  eurfhoe,  rnd  D^  -  x  «ino»  the  nediuw  t«  homogeneous 
the  diffracted  ray  io  a  eurfhoe  geodesic,  lies  along  the  unit  oornsl  to  the 
ray,  and  lias  along  the  biannual*  At  a  point  P  on  the  straight  diffracted 
ray  Milch  leaves  the  surface  at  ««  dafius  the  triad  of  vectors 

•  The  following  discussion  la  adapted  froa  cacti  on  f*  or  fen.  The  analogous 
theory  for  an  lnbonogeneous  aediua  baa  nob  yet  bean  devalued. 
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by  ssttiac  DV(P)  -  Dy(P^).  We  denote  the  noapononts  of  the  electric  neli 
vector  £  **  eiks2  in  the  three  directions  Dy  hy  £  y. 

We  new  assume  that  g,  ,  is  zero  and  that  the  components  and 
g,  propagate  independently  of  each  other  end  satisfy  the  equations  of 
Chap '.nr  A.  Then  from  (Aie-ij) 


Here  £  *(0^)  -  S>1(Q1),Dy(Q1)  is  the  cc^cnent  of  the  Incident  field  £*  at  * 

^  in  the  direction  Dw,  end  the  other  (juentities  ere  defined  in  Section  Alfl.  ’ 

a 

The  diffraction  coefficients  djy  end  the  decay  exponents  o1y  ere  different 

for  the  two  oospeneota.  At  P,  the  difiVuatvd  field  associated  with  e  diffracted  • 

ray  ie  given  by  (tS)  and 

i 

s"(M .  ijtnojtn  *  s!(f)i>3(»).  (k i) 

nr  .  prrreetlir  ms-tnt  -oth  »)«t  tb.  tj}  «t 

Ojj  for  the  tengrntlal  cccymant  £  ^  are  the  ten*  ea  vhoee  for  a  eoelsr  field 
aetiafyinq  the  condition  u  «  o  on  the  surface  and  henoo  can  be  obtained  fns 
(A19«13  .  11)  by  setting  «  »•.  Slailarly  the  coefficients  for  the 

sonel  corpouent  £  ?  are  the  sane  as  those  tor  a  sealer  field  mtlafylne  *-U* 
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roruUtior  r-  *  0  or.  the  surface  and  hence  cun  be  obtained  ftroa  the  suae 
o* 

equations  by  setting  z  --  0.  0;'  course  fa  mist  also  cat  *(X)  •  0  and  n(X)  •  n  in 
then*  equations  since  the  medium  is  homogeneous. 
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